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ABSTRACT 



Context. The present number of Galactic Open Clusters that have high resolution abundance determinations, not only of [Fe/H], but 
also of other key elements, is largely insufficient to enable a clear modeling of the Galactic Disk chemical evolution. 
Aims. To increase the number of Galactic Open Clusters with high quality measurements. 

Methods. We obtained high resolution (R~30000), high quality (S/N~50-100 per pixel), echelle spectra with the fiber spectrograph 
FOCES, at Calar Alto, Spain, for three red clump stars in each of five Open Clusters. We used the classical Equivalent Width analysis 
method to obtain accurate abundances of sixteen elements: Al, Ba, Ca, Co, Cr, Fe, La, Mg, Na, Nd, Ni, Sc, Si, Ti, V, Y. We also 
derived the oxygen abundance through spectral synthesis of the 6300 A forbidden line. 

Results. Three of the clusters were never studied previously with high resolution spectroscopy: we found [Fe/H]=+0.03±0.02 
(±0.10) dex for Cr 110; [Fe/H]=+0.01±0.05 (±0.10) dex for NGC 2099 (M 37) and [Fe/H]=~0.05±0.03 (±0.10) dex for NGC 2420. 
This last finding is higher than typical recent literature estimates by 0.2-0.3 dex approximately and in better agreement with Galactic 
trends. For the remaining clusters, we find: [Fe/H]=+0.05±0.02 (±0.10) dex for M 67 and [Fe/H]=+0.04±0.07 (±0.10) dex for 
NGC 7789 . Accurate (to ~0.5 km s"') radial velocities were measured for all targets, and we provide the first high resolution based 
velocity estimate for Cr 1 10, <V,.>=41.0±3.8 km s"' . 

Conclusions. With our analysis of the new clusters Cr 110, NGC 2099 and NGC 2420, we increase the sample of clusters with high 
resolution based abundances by 5%. All our programme stars show abundance patterns which are typical of open clusters, very close 
to solar with few exceptions. This is true for all the iron-peak and s-process elements considered, and no significant or-enhancement 
is found. Also, no significant sign of (anti-)correlations for Na, Al, Mg and O abundances is found. If anticorrelations are present, 
the involved spreads must be <0.2 dex. We then compile high resolution data of 57 OC from the literature and we find a gradient 
of [Fe/H] with Galactocentric Radius of -0.06±0.02 dex kpc ', in agreement with past work and with Cepheids and B stars in the 
same range. A change of slope is seen outside i?GC=12 kpc and [ff/Fe] shows a tendency of increasing with Rqq. We also confirm the 
absence of a significant Age-Metallicity relation, finding slopes of -2.6±I.I 10"" dex Gyr"' and I.I±5.0 10"" dex Gyr"' for [Fe/H] 
and [ff/Fe] respectively. 
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^ 1. Introduction 

Open clusters (hereafter OC) are the ideal test particles in the 
study of the Galactic disk, providing chemical and kinemat- 
ical information in different locations and at different times. 
Compared to field stars, they have the obvious advantage of 
being coeval groups of stars, at the same distance and with a 
homogeneous composition. Therefore, their properties can be 
determined with smaller uncertainties. Several attempts have 
been done in the past to derive two fundamental relations us- 
ing OC: the metallicity gradient along the disk and th e age 



metal l icity relation (hereafte r AMR) of th e disk (e.g. iJanes , 



1979; 



2002t 



Panada & To"sil 119801; iTwarog etalL 119971; iFriel et al . 
Chen et all 120031: ISalaris et all |2004 . but they were 



Send offprint requests to: E. Pancino 

* Based on data collected with the fiber spectrograph FOCES at the 
2.2m Calar Alto Telescope. Also based on data from the 2MASS survey 
and the WEBDA, VALD, NIST and GEISA online databases. 



hampered by the lack of large and homogeneous high quality 
datasets. 

In particular, the lack of a metallicity scale extending to so- 
lar metallicity w ith comparable prec i sion to the lower metallic - 
ity regimes (i.e., IZinn& Wesi Il984l: ICarretta & GrattonL fT997h 
represents the main problem from the point of view of (i) the 
study of the Galactic disk; ( ii) tests of stellar evolution models 
for younger and more metal-rich simple stellar populations and 
( Hi) the use of those stellar populations as templates for extra- 
gala ctic studies of pop ulation synthesis. Of the ~1700 known 
OC dPias et all I2OO2I and updates), only a subset of -140, 
i.e., 8% of the total, possesses some metallicity determination. 
IVlost of these have been obtained through different photomet - 
ric studies in the Washington (e.g. i Geisler et al .l '1991', 1 19921). 
DDO (e.g. [C laria et al., 1999), Stron igren (e. g. Brunt t et all 
11999'; Twarog et al., 2003), UBV (e.g. ICameronl il985) and IR 
(e.g. Tiede et al., 1997) photometric systems and passbands, giv- 
ing often rise to consi derable differen ces with those obtained 
from spectroscopy (see iGrattonl l200(i and references therein). 
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Table 1. Observing Logs and Programme Stars Information. 



Cluster 


Star 


Q'J2000 

(hrs) 


^72000 

(deg) 


B 

(mag) 


V 

(mag) 


Ic 

(mag) 


R 

(mag) 


K5 
(mag) 


^exp 


Mar) 

(sec) 


S/N(^6000A) 


Cr 110^ 


2108 


06:38:52.5 


+02:01:58.4 


14.79 


13.35 


— 


— 


9.76 


6 


16200 


70 




2129 


06:38:41.1 


+02:01:05.5 


15.00 


13.66 


12.17 


12.94 


10.29 


7 


18900 


70 




3144 


06:38:30.3 


+02:03:03.0 


14.80 


13.49 


12.04 


12.72 


10.19 


6 


16195 


65 


NGC 2099 (M 37J^ 


067 


05:52:16.6 


+32:34:45.6 


12.38 


11.12 


9.87 


— 


8.17 


3 


3600 


95 


148 


05:52:08.1 


+32:30:33.1 


12.36 


11.09 


— 


— 


8.05 


3 


3600 


105 




508 


05:52:33.2 


+32:27:43.5 


12.24 


10.98 


— 


— 


7.92 


3 


3900 


85 


NGC 2420^ 


041 


07:38:06.2 


+21:36:54.7 


13.75 


12.67 


11.61 


12.13 


10.13 


5 


9000 


70 




076 


07:38:15.5 


+21:38:01.8 


13.65 


12.66 


11.65 


12.14 


10.31 


5 


9000 


75 




174 


07:38:26.9 


+21:38:24.8 


13.41 


12.40 






9.98 


5 


9000 


60 


NGC 2682 (M 67^ 


0141 


08:51:22.8 


+ 11:48:01.7 


11.59 


10.48 


9.40 


9.92 


7.92 


3 


2700 


85 


0223 


08:51:43.9 


+ 11:56:42.3 


11.68 


10.58 


9.50 


10.02 


8.00 


3 


2700 


85 




0286 


08:52:18.6 


+ 11:44:26.3 


11.53 


10.47 


9.43 


9.93 


7.92 


3 


2700 


105 


NGC 7789^ 


5237 


23:56:50.6 


+56:49:20.9 


13.92 


12.81 


11.52 




9.89 


5 


9000 


70 




7840 


23:57:19.3 


+56:40:51.5 


14.03 


12.82 


11.49 




9.83 


6 


9000 


75 




8556 


23:57:27.6 


+56:45:39.2 


14.18 


12.97 


11.65 




10.03 


3 


5400 


45 



" Star names and Ic & R magnitudes from|D awson & lannal ( Il998h : coordinates and B & V magnitudes from lBragaglia & Tosil ( l2003h : 
K5 magnitudes from 2MASS 

'' Star names from 'van Zeipel & Lindgre3 ( 1192 Ih : coordinates from lKiss et al.l ( 1200 Ih : B (& V magnitudes from lKalirai et al.l ( 120011) : Ic 
magnitudes from Nilakshi & Sagar ( 2002); Ks magnitudes from 2MASS 

Star names from Cannon & Lloyd (1970); coordinates from Stetson (2000) and i Lasker et alj dl990l) ; B & V magnitudes from 
lAnthonv-Twarog et al. ( 1990); Ic & R magnitudes from Stetson (2000); magnitudes from 2MASS 

Star names from Fagerholm ( 1906); coordinates from iFan et alJ ([1996); B, V c& Ic magnitudes from lSandquisd ( l2004h ; B magnitude for 
star 286 and R magnitudes from janes & Smithltl984i) ; K5 magnitudes from 2MASS 

Star names and V & L magnitu des from lGim et alj U998a); J1950 coordinates from lKustnej il923h . precessed to 12000; B magnitudes 
from lMocheiska & Kaluznvl ( 119991) ; K5 magnitudes from 2MASS 



In a much smaller number of clusters, abundan ces have been 
derived from low-resolution spectroscopy (e.g., ICarrera et al.L 
120071; IWarren & CoR 12008 ). with admir able attempts to ob- 
tain large and ho mogeneous datasets (see iFriel & Janesl 119931; 
iFriel et alll2002h . in spite of the non-negligible uncertainties in- 
volved in the procedure. 



A few research groups (see Section |6] for more details) are 
presently obtaining high quality spectra and are producing more 
precise abundance determinations. The study of elements other 
than the iron-peak ones (such as a, s- and r-process, light ele- 
ments), allows one to put more constraints on the sites of produc- 
tion of those elements (SNe la, SNe II, giants and supergiants, 
Wolf-Rayet stars), and therefore on their production timescales. 
These are fundamental ingredien ts for the chemical evolution 
modeling of the Galactic Disk dTosil 119821 ; IChiappini et al.L 
l200lHColavitti et al.Ll2009l) . 



For these reasons, we obtained high resolution spectra for a 
sample of poorly studied old OC. We present here the detailed 
abundance analysis of five clusters observed during our first run 
at Calar Alto. Observations and data reductions are described in 
Section |2l the linelist and equivalent width measurements are 
detailed in Section[3]while the abundance analysis methods and 
results are presented in SectionlU abundance results are then dis- 
cussed and compared with literature results in Sections |5j|6] and 
|7j finally, we summarize our results and draw our conclusions in 
SectionlH 



2. Observational Material 

Three red clump starfl were selected in each of the tar- 
get clusters using the WEBDAQ database dMermilhodl [19951) 
and the 2MASSQ survey data for the infrared Ks magnitudes 
(Skrutskie et al., 2006). More details on references for star 
names, coordinates and magnitudes can be found in Table [1] 
while the position of our targets in the Color Magnitude 
Diagrams (CMDs) obtained from WEBDA are shown in 
Figure [1] 

Observations were done between the 1st and 10th of January 
2004 with the fiber echelle spectrograph FOCES at the 2.2 m 
Calar Alto Telescope, in Spain. The sky was generally clear, 
though a few nights had thin cirrus and sometimes clouds, forc- 
ing us to increase the exposure times considerably. All stars were 
observed in 3-16 exposures lasting 15-90 min each, depending 
on the magnitude, until a global S/N ratio between 70 and 100 
(per pixel) was reached around 6000 A (Table[TJ. Each night we 
took a sky exposure lasting as the longest one of the night. The 
sky level was negligible for all exposures with S/N>20, therefore 
exposures with S/N<20 were neglected and we did not subtract 
the sky, to avoid adding noise to the spectra. Sky emission lines, 
even in the red part of the spectrum, did only occasionally disturb 



' Though fainter than the brigthest giants, red clump stars have the 
advantage of a higher gravity and temperature, that reduces consider- 
ably line crowding. Also, clump stars are easy to identify even in the 
sparsest cluster, maximizing the chance of choosing cluster members. 

- http://www.univie.ac.at/webda 

^ http://www.ipac.caltech.edu/2mass 2MASS (Two Micron 
All Sky Survey) is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 
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Fig. 1. V, (B-V) Color Magnitude Diagrams of the programme clusters (from the WEBDA), with the location of our target stars. 



the measurements of some absorption Unes, that were discarded. 
The spectral resolution was R ^ 30 000 for all spectra. 

2.1. Data Reductions 

Data reductions were done with within the imred and 

echelle packages. The following steps were performed: bias sub- 
traction, flatfielding, order tracing with optimal extraction, back- 
ground subtraction, wavelength calibration with the help of a 
Thorium- Argon lamp, and final merging (and rebinning) of over- 
lapping orders. The one-dimensional spectra obtained from dif- 
ferent exposures (with S/N>20) were median-averaged to pro- 
duce one single high S/N spectrum for each star, used for equiv- 
alent width measurements (Section [33] l. Finally, the noisy ends 
of each combined spectrum were cut, allowing for an effective 
wavelength coverage from 5000 to 9000 A. 

Sky absorption lines (telluric bands of O2 and H2O) were 
removed using the IRAF task telluric with the help of two hot, 
rapidly rotating stars, HR 3982 and HR 8762, ch osen from the 
Bright Star Catalogue (Ho ffleit & Jascheklll99l1) . HR 3982 and 
HR 8762 were observed each night at an airmass not too dif- 
ferent from the scientific targets. Residuals of the correction in 
the red part of the spectrum (for example from the strong O2 
band around 7600 A) prevented us from using most of the corre- 
sponding spectral regions for our abundance analysis. Also, af- 
ter 8400 A, the echelle orders do not overlap anymore and small 
gaps appear. 

2.2. Radial Velocities 

Radial velocities were m easured with the help of DAOSPEC 
(IStetson & Pancinol 1200 8'. see also Section IX3T l. Measurements 
were based on ^360 absorption lines of different elements (see 
SectionO with typical measurement errors on the mean of about 



Image Reduction and Analysis Facility. IRAF is distributed by the 
National Optical Astronomy Observatories, which is operated by the 
association of Universities for Research in Astronomy, Inc., under con- 
tract with the National Science Foundation 



Table 2. Heliocentric radial velocities measurements and Icr er- 
rors (Vr + SVr)here for cach programme star. Literature measure- 
ments, when available, are also reported with their uncertainties 

(Vr + 6Vr)m. 



Cluster 


Star 


(Vr±SVr)l„r, 


(V,- ± crVrh, 






(km s-') 


(km s-') 



Cr 110 


2108 


44.47 = 


tO.54 


45 H 


t 8^ 




2129 


38.74 = 


tO.64 


45 H 


t 8~ 




3144 


37.81 = 


tO.69 


45 H 


t 8° 


NGC 2099 (M 37) 


067 


8.79 ± 


0.10 


8.04 ± 


0.19^ 




148 


9.05 ± 


0.36 


8.73 ± 


0.19~ 




508 


9.40 ± 


0.28 


8.96 ± 


0.19* 


NGC 2420 


041 


74.23 = 


tO.87 


79.0 = 


t4.2_^ 




076 


75.49 = 


t041 


75.7 = 


t 8.6~ 




174 


73.66 = 


t 1. 17 


68.0 = 


tO.9'- 


NGC 2682 (M 67) 


141 


35.23 = 


t0.05 


33.5 = 


t2.3^ 




223 


34.92 = 


t0.31 


33.0 H 


t I.l" 




286 


38.90 = 


tO.46 


33.2 H 


t 2.0'' 


NGC 7789 


5237 


-51.06 


±0.99 


-57.17 


± 0.23^ 




7840 


-49.06 


± 0.81 


-49.21 


± 0.28~ 




8556 


-53.37 


±0.81 


-54.10 


± 0.32" 



" Cluster average bv lCarrera et based on 8 stars. 

* iMermilliod et alJ j2008h . su persed ing Mermilliod et al](ll996l). 
'^ Aver age of measu rements bv lLiu & JanesI ( ll987l) . IFriel et al.l il989l) 
and iScottetalJjl995h . 

'' Average of mea s urements by iMathieu et aP 
Pilachowski et al.' (1 19881) . iFriel et al.l H9m. iFriel & JaM | 



1993), 



Suntzeff et al. ( 1992), Sco tt et aLUl'99a) . lMelo et al. (2001:). lYong et all 
j2005). Yadav et al. (2008). 
' FOim et al. U998b,) . 



0.1 km s ^ All measurements were performed separately on the 
one-dimensional spectra extracted from the single exposures for 
each star, including those with S/N<20, that were not used for 
the abundance analysis. In this way, we could check that no sig- 
nificant radial velocity variations were present. 
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Table 3. Stellar Parameters for the programme stars. 



Cluster 


Star 


(K) 


(K) 


log g^P''°'> 

(cgs) 


logg('P<-c) 
(cgs) 


(km s-') 


(km s-') 


(Mo) 


Cr 110 


2108 


4914±111 


4650 


2.32±0.16 


2.8 


1.20±0.02/1.62±0.04 


1.4 


1.9+0.1 




2129 


5056±223 


4950 


2.53±0.19 


2.7 


1.17±0.02/1.52±0.06 


1.4 


1.9+0.1 




3144 


5112±258 


4800 


2.49±0.20 


2.8 


1.18+0.02/1.49+0.08 


1.3 


1.9+0.1 


NGC 2099 (M 37) 


067 


4773±119 


4550 


2.15±0.32 


2.7 


1.22+0.03/1. 68+0.05 


1.5 


2.7+0.2 




148 


4708±116 


4550 


2.10±0.31 


2.7 


1.28+0.03/1.72+0.04 


1.5 


2.7+0.2 




508 


4715±115 


4500 


2.05±0.32 


2.8 


1.23+0.03/1.72+0.04 


1.5 


2.7+0.2 


NGC 2420 


041 


4616±77 


4850 


2.43±0.06 


2.6 


1.19+0.01/1.76+0.03 


1.4 


1.4+0.1 




076 


4755±103 


4800 


2.51±0.06 


2.6 


1.18+0.01/1.69+0.03 


1.6 


1.4+0.1 




174 


4730±77 


4800 


2.39±0.05 


2.6 


1.19+0.01/1.71+0.03 


1.5 


1.4+0.1 


NGC 2682 (M 67) 


141 


4590±100 


4650 


2.42±0.04 


2.8 


1.19+0.01/1.78+0.04 


1.3 


1.3+0.1 




223 


4594±100 


4800 


2.46±0.04 


2.8 


1.18+0.00/1.78+0.04 


1.3 


1.3+0.1 




286 


4653±103 


4850 


2.45±0.04 


2.8 


1.18+0.01/1.75+0.04 


1.4 


1.3+0.1 


NGC 7789 


5237 


4868±168 


4900 


2.53±0.15 


2.8 


1.17+0.02/1.63+0.04 


1.2 


1.8+0.1 




7840 


4759±131 


4800 


2.47±0.15 


2.7 


1.18+0.02/1.69+0.04 


1.5 


1.8+0.1 




8556 


4775±135 


4900 


2.54±0.15 


2.9 


1.17+0.02/1.68+0.04 


1.4 


1.8+0.1 



Heliocentric corrections were computed with the IRAF 
task rvcor, which bears a negligible uncertainty of less than 
0.005 km s"' . Since we did not observe any radial velocity stan- 
dard and our calibration lamps were not taken simultaneously, 
we used telluric absorption lines to find the absolute zero point 
of our radial velocity measurements. In particular, laboratory 
wavelengths of the H2O absorption bands around 5800, 6500, 
7000, 7200, 8000 and 8900 A and of the O2 absorption bands 
around 6 300, 6900 and 7600 A were obtained from the GE1SA0 
database (iJacquinet-Husson et al.L[T999ll2005h and we measured 
their radial velocity in our programme stars. The resulting zero- 
point corrections, based on 200-250 telluric lines, amount gen- 
erally to no more than ±1 km s ' with typical errors on the mean 
of about 0.5 km s"', approximately five times larger than those 
on the radial velocity measurements. 



After applying the above corrections, and propagating the 
corresponding uncertainties, we computed a weighted average 
of the heliocentric velocities estimates for each exposure (see 
Table |2]i. All the programme stars appear to be radial velocity 
members of the observed clusters, with the possible exception 
of star 2108 in Cr 110, that has a slightly higher velocity than 
2129 and 3144. However, since the value for 2108 is within 3cr 
from the mean value for the cluster, we decided not to reject this 
star We can provide the first radial velocity estimate based on 
high resolution for Cr 110: <V',.>=41.0±3.8 km s Our deter- 
minations are generally in good agreement with literature val- 
ues within 3cr, except maybe for star 5237 i n NGC 7789, which 
is marginally dicrepant with the estimate bv lGim et al.l (Il998bl) . 
However, there is perfect agreement between the two studies for 
the other two stars of NGC 7789, and our estimate appears more 
in Une with a membership of 5237. In conclusion, we consid- 
ered all the programme stars as likely radial velocity members 
of their respective clusters. 



2.3. Photometric Parameters 

We first computed the dereddened color^ (B-V)o, (V-Ic)(JZl 
and (V-K7-C5)(0. The adopted E(B-V) values are indicated 
in Table |4] wh ere ECV-Ic ) was obtained with the red- 
dening laws bv iDeanetal'l (1 19781) . and E(V-K7-C5) with 
ICardelh et all (Il989h . We were then able to obtain Tetf and 
the bolometri c correction BCy for each programme star, us- 
ing both the 'Alonso et al.' ('1999) and the (theoretical and 
empyrical) Montegriffo et al. (1998) color-temperature rela- 
tions, taking into account the uncertainties on magnitudes 
and reddening estim ates. Th e average difference bet ween the 
lAlonso et all d 19991) and the iMontegriff'o et alj ( Il998h temper- 
atures w as Areff=178± 66 K (for the empyrical calibration of 
Montegri ffo et al.Lll998h and Areff=127±6 9 K (for the theoreti- 



Montegriffo et~ail fT998h . We averaged all the 



cal calibration of 
above Teft estimates to obtain our photometric reference values 
with their Icr uncertainties (TableO. 

Gravities were obtained from Tej and BCy using the funda- 
mental relations 



logA.iog^.21og| 



0.4(Mboi - MboLo) 



-4 log -A +2 log ^ 

eff,0 " 



' http : // ara . Imd . polyte chnique . fr/htdOC-public/pro | 

ducts/GEISA/HTML-GEISA/ 



where red clump mas ses were derived using Table 1 of 
iGirardi & SalarisI (1200 ll) . and are also shown in Table [3] We as- 
sumed log go = 4.437, Teff.o = 5770 K and Mboi.o = 4.75, i n 
conformity with the lAU recommend ations ( Andersenl 1 19991) . 
The difference between the lAlonso et a l. ( 1999) and t he (em pyri- 
cal and theoretical, respectively) [Mon tegriffo e t al.l (1 19981) esti- 
mates was Alogg=0.20+0.13 and Alogg=0.18±0.13. As above, 
we averaged all our estimats to obtain photometric gravities 
log g(pi""'i (TableO and their Icr uncertainties. 

Finally, a photometric estimate of the microturbulent ve- 
locities V, was obtained using both the prescriptions of 

* Since R magnitudes are available for less than half of our sample, 
we decided to ignore them in the following. 

' After dereddening, (V-Ic) was also converted into (V-I/) using the 
prescription by Bessell ( 1979), to be used with the color-temperature 
relations by Alonso et al. ( 1999.) . 

** We computed the KrCS magnitudes from the 2MAS S TCS magni- 
tudes using the prescription bv lKirmfiari &Castellil(l2002h . 
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Table 4. Input cluster parameters. 



Cluster 


E(B-V) 


(m-M)v 


log Age 




(mag) 


(mag) 


(dex) 


Cr 110^ 


0.54±0.03 


13.37±0.38 


9.23±0.15 


NGC 2099 (M 37]^ 


0.27±0.04 


11.53±0.19 


8.61±0.16 


NGC 2420;. 


0.04±0.03 


11.88±0.27 


9.47±0.17 


NGC 2682~(M 67^;;^ 


0.04±0.02 


9.67±0.11 


9.64±0.05 


NGC 7789'' 


0.27±0.04 


12.23±0.20 


9.21±0.12 



" Averages of measurements by Tsarevskii & Abakumovl il97lh . 
iDawson & Imri3 ( (1998) and,Braga glia & Tosi (200^). 

Averages o f mea s urements bv | Westl (fT967ah. ITerm ens & Heifer 
(Il975h. iLynpl (Il987 l). 'Mermilliod et al.' |'199^. 'Twarog et al. (1997), 

(2001), Nilakshi & Sagar (2002), 
iBragaglia & Tosii (|200Q) and 



Kiss et al.l T2001I). Kalirai et al 
Grocholski & Saraiedinii (I2003t) , 
Hartmanetal. (200^~^ 

Averages of meas ur ements by 
Barbaro & Pigattg| (|1984|), |Vandenber gl 



McClure et alj 
19851), iLyngal 



Anthony-Twarog et al. 
Demarque et al. ( 1994^ 



( ^1974) . 
09871), 
(1994), 



^4i)..iwarog 

Grocholski & Sarai edinj d2003l) . 

Averages of measurements by 'Racing 



(1990), Carraro & Chiosi 
Twarog et al.i 11997.) . . Tadross et aL (,200Z) and 



19711) . IB arbaro & Pigatto! 
1984|), IVandenberg (1985), Nissen et air il987b. [Hobbs & Thorburn 
19911), Gilliland & Brown (1992), Demarque et al. (1992), 
i Meynet et al. (1993), Montgomery et al. (1993), Carraro & Chiosi 
( 1994), Dinescu et al. ( 1995)^ Fan et al. ( 1996), Twarog et al. ( 1997), 
[Saraiedin i (1999), Tadross et al. C2002). ^Grocholski & Saraiedini 
(120031), IV andenBe rg & ClemI ( l20"03h . iLaugalvs et alj d2004h and 
ISandguislI ( 12004.) . 

Averag e s of measurements b y iBurbidge & Sandagd dl958t), 

lArcl (I196Z ). 'Strom & Strom' ("1970), 'Jennens & Heifer' ("1975^, 'Janes' 
19770, IC laria (1979), Twarog c& Tyson (1985), Lynga (1987), 



|m azzei & Pigatto (1988), Friel & Janes (1993), Martinez Roger et al, 
( 119941) . 'Janes & Phelps (1994), Twarog et al. (199 7), ,Gim et_aL 
( 1998b), Saraiedini ( 1999), Vallenar i et al.li'lOOOh . iTadross et aljj2002l) ' 
and lBartasiute & Tautvaisiene. ( .20041) . 



Ramfrez & CohenI (l2003l) . v, = 4.08 - 5.01 10 Teif, and of 
Carretta et all (1200 4'), v, 1.5-0.13 log^. The latter takes into 



account the systematic effect discussed by Magain (1984) and 
is on average lower by Avf - 0.49 ± 0.08 km s ' than the one 
bv lRamfrez & CohenI (12003]). Ho wever, the actual amount of the 
correction for the MagainI ( 1984| ) effect depends heavily on the 
data quality (i.e., resolution, S/N ratio, number of lines used, 
log gf values etc.). Therefore we chose not to average the two 
estimates, but to use them as an indication of the (wide) v, range 
to explore in our abundance analysis (see Section l^TI ). 



3. Linelist, Atomic Data and Equivalent Widths 

We created a masterlist of absorption lines by visually compar- 
ing our spectra with the the UVES solar spectrurrQ in the range 
5000-900 A, and with the lin e lists extracted fron i the VALC[lB 
datab ase dKupka et all Il999h and the MooreQ dMoore et all 
11966 1) solar atlas. The masterlist was fed to DAOSPEC and EW 
were measured for all our programme stars. A first selection was 
applied to reject all those lines that were measured in 10 stars 
or less (out of 15) and that had EW systematically larger than 
250 mA. Later, after performing a first crude abundance analy- 
sis (see Section|4]T), we rejected all the lines that gave sistemati- 



' http : //www ■ eso ■ org/ obser ving/df o/quality/UVES/pipel | 
ine/solar_spectrum . html 

jhttp://www. astro. uu.se/~vald7] 
" |f tp : //ftp . noao . edu/f ts/linelist/Moore | 



cally discrepant abundances, especially if the formal DAOSPEC 
relative error (5EW/EW, Figure|2]i was around 15% or more, and 
the DAOSPEC quality parameter was above 1 .5 (for more details 
about DAOSPEC error and quality parameter, see Section [33T l. 
The final linelist, including atomic data and EW measurements 
for all programme stars, contains 358 absorption lines of 17 
species, and can be found in the electronic version of Table |5] 
Atomic data include laboratory wavelengths, excitation poten- 
tials and log gf values, which are always taken from the VALD 
database, with the exceptions listed below. 

3.1. a Elements Atomic Data 

The only (T-element for which we had clear problems with the 
atomic data was magnesium. The lines with;t'„=5.75 eV (7060 
and 7193A) gave discrepant abundances by ~1.5 dex with re- 
spect to the average of all Mg lines. We compared our VALD 
log gf with the NIST0 database of atomic data and noticed a 
difference of 1.4 dex for the Xex-^.IS eV lines, while all the 
other Mg I had very similar log gf values in both databases. The 
NIST log gf values abundances of the Xex-5.15 eV lines gave 
abundances in much better agreement with the other Mg I lines 
and with the literature Mg abundances for OC, therefore we used 
the NIST values for those lines, instead of the VALD ones. 

Another element with uncertain log gf values is Calcium. 
As an example, for the 9 lines we use, there is an average 
difference of logg/Af/s7-log^/vALD=-0.17+0.18 dex, which is 
not statistically significant given the large cr. Also, the NIST 
\oggf values for those 9 lines all range from D to E, which 
means that they are largely uncertain. Finally, our solar abun- 
dance (Sectiong^ll gives [Ca/H]=-0.09+0.03 (±0.03) dex if we 
use the VALD loggf and [Ca/H]=-H0.08±0.03 (±0.03) dex with 
the NIST ones, which is equally compatible with zero within 3 
sigma. Summarizing, there is large uncertainty on the Calcium 
loggf determinations, and we really should keep in mind that 
there is an additional ~0.2 dex uncertainty on all [Ca/Fe] deter- 
minations in the literature. 

For the synthesis of the [O I]-Ni I blend at 6300A, we 
used the VALD log gf for oxygen, but we chose to use the 
[Johansson etal] (l2003h log gf for Ni I at 6300.35A, which is 
lower (-2. 1 1 dex instead of -1 .74) and gives oxygen abundances 
more in line with the other a-elements. 



3.2. Heavy Elements Atomic Data 

For neodymium, we could only find three reliable lines, which 
apparently do not need any detailed HFS (hyper-fine splitting) 
analysis (i Aoki et al.Ll200l . at 5092, 5249 and 5485 A. However, 
the spread of their abundance s was quite high (T able[7li. The lab- 
oratory log gf published by iDen Hartog et al.l (i2003) are very 
similar to the ones from VALD, except for the 5485 A line, 
where they di ffer by .14 dex. Therefore, since the log gf values 
by [ben Hart og et al.l (12003) slightly reduced the spread in the 
[Nd/Fe], we used them instead of the VALD ones (see Table|5]l. 



3.3. Equivalent Widths with DAOSPEC 

The full description of how DAOSPEC works, including com- 
parisons with the literature and several ex periments with ar- 
tificial and real spectra, can be found in Stetso n & Pancinol 
(120081) . The instructions on how to install, configure and use 



jhttp : //physics . ni st . gov/PhysRefPata/ ASD/index . html j 
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Table 5. Equivalent Widths and Atomic Data of the programme stars. The complete version of the Table is available at the CDS. 
Here we show a few lines for guidance about its contents. 



Cluster 


Star 


A 


Elem 


/Vex 


logg/ 


EW 


SEW 


Q 






(A) 




(eV) 


(dex) 


(mA) 


(mA) 




Cr 110 


2108 


5055.99 


Fel 


4.31 


-2.01 


41.2 


3.3 


0.371 


Cr 110 


2108 


5178.80 


Fel 


4.39 


-1.84 


45.4 


7.1 


0.851 


Cr 110 


2108 


5294.55 


Fel 


3.64 


-2.86 








Cr 110 


2108 


5285.13 


Fel 


4.43 


-1.64 


50.1 


5.2 


0.607 


Cr 110 


2108 


5295.31 


Fel 


4.42 


-1.69 


38.3 


9.5 


1.958 
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Fig. 2. The behaviour of DAOSPEC relative errors 5EW/EW 
versus EW is shown (top panel) and the 15% error limit is 
marked with a dotted line. The behaviour of the quality parame- 
ter Q is shown versus EW (middle panel) and wavelength (bot- 
tom panel) and the Q=1.5 limit is marked with a dotted line. See 
Section [33] for more details. In all panels, grey dots represent 
all the lines measured by DAOSPEC, while black dots represent 
lines cross-identified with our input linelist. 



ter line list. Since our spectra are rebinne linear in wavelength, 
we scaled the FWHM with A. Figure |2] shows their behaviour 
dEW is the formal error of the gaussian fit that DAOSPEC out- 
puts, and 6EW/EW can be used to select good measurements, 
since smaller lines are noisier and tend to have higher relative 
errors. The quality parameter Q, instead, is the result of the com- 
parison of local residuals around each line with average residuals 
on the whole spectrum. As a result Q tends to be worse for strong 
lines, because the Gaussian approximation does not hold so well 
anymore. Also, Q gets worse at the blue side of the spectrum, 
where the S/N ratio is lower. In the region around 7700 A, where 
the residuals of the prominent O2 telluric band disturb the mea- 
surements, Q reaches its maximum. The measured EW for our 
programme stars are shown in the electronic version of Table |5] 
along with the 6EW and Q parameter estimated by DAOSPEC. 



3.4. EW uncertainties 

We used the formal errors of the Gaussian fit computed by 
DAOSPEC only to reject bad measurements from our initial 
line list. The actual abundance errors due to the EW measure- 
ments process itself were instead computed later, as explained in 
Section |43] 

To compute the EW uncertaint y due to the continuum place- 
ment, we used Equation 7 from IStetson & Pancino (2008) to 
derive the effective uncertainty on the continuum placement 
(±AC/C), which turned out to be significantly smaller than 1%. 
We first lowered the "best" continuum by AC/C and measured 
EWs again, obtaining EW(_), then we increased it by the same 
amount and measured EW(+). The differences with the "correct" 
EW measurements, AEW(_) and AEW(+) were averaged to pro- 
duce the actual AEW for each line. The typical resulting uncer- 
tainty, due only to the continuum placement, was approximately 
constant with EW and of AEW ^l mA approximately (see also 
Figure 2 by Stetson & Pancind.12008) . Such a small uncertainty 
was neglected because it had a much smaller impact on the re- 
sulting abundances than other sources of uncertainty considered 
in Sections l43]andl4!4l 



DAOSPEC can be found in "Cooking with DAOSPEC'^ In 
short, DAOSPEC is a Fortran program that automatically finds 
absorption lines in a stellar spectrum, fits the continuum, mea- 
sures EW, and identifies lines with a laboratory linelist; it also 
provides a radial velocity estimate (Section |2!2] ). 

As described in Section[3] we used the DAOSPEC errors and 
quality parameter, Q, to select good absorption lines in our mas- 



http : //www3 ■ cadc -ccda . hia-iha .nrc-cnrc . gc . ca/commun 
ity/STETSON/daospec/ ; |http://www.bo.astro.it/~pancino 
/pro jects/daospec. html 



3.5. Comparison witli Literature EW 

To our knowledge, on ly o ne of our target stars was studied before 
by lYongetaLl(l2005h and lTautvaisiend (l2000h . with a resolution 
and S/ N similar to ours, i.e. star 141 in M 67. While l^ng et al.l 
(2005) do not publ ish their E W measurements, we can compare 
with the ones bv Tautvaisiene (2000). The authors provided two 
sets of EW, the former derived from a spectrum with R^30000, 
and the latter from a spectrum with R^60000. We have 48 lines 
in common with the R^30000 set and 36 with the R^60000 set. 

Figure [3] shows good agreement between our EWs and the 
R-30000 set. We found just a sistematic offset of EWroo- 
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Table 6. Literature comparison for star 141 in M 67. 
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Fig. 3. Comparison of our EW measurements with iTautvaisiend 
for star 141 in M 67. The top panel shows the compar- 
ison of DAOSPEC EW with the R^30000 set by ITautvaisiend 
(120001 based on 48 lines in common. The middl e panel show 
the sam e comparison, but for the R^60000 set by ITautvaisiend 
(1200 0*). based on 36 lines. The bottom panel shows the 15 lines 
in common between the R^BOOOO and the R^60000 EW sets, 
considering only the lines included in the two upper panels. 



DAOSPEC - 4.0 ±7.3 mA, which corresponds to a contin- 
uum placement difference of about 1% (see also Section U4b . 
A possible trend with EW was visible when comparing our 
EWs with the R^60000 set, with no systematic offset (AfiW - 
-0.3 ± 9.3 mA). On the one hand, this means that our continuum 
placement ag rees much better with the R^60000 continuum by 
lTautvaisiener (i2000) than with the R^30000 one. On the other 
hand, we noted that a possib le trend is also visible when com- 
paring the ITautvaisiend ( 1200 0) measurements at R^30000 with 
those at R^60000. In conclusion, we considered our EW mea- 
surements in good agreement with the Tautvai siend (l2000l) ones, 
given the involved uncertainties (see also Table|6]l. 



4. Abundance Analysis 

4.1. Best Model Search 

A preliminary abundance determination was done using the 
photometric parameters (Section |4]i, which allowed us to iden- 
tify and remove those lines in our list that gave systematically 
discrepant abundances. We found largely discrepant Fe I and 
Fe II values when using the photometric parameters (by ^0.5- 
0.9 dex), indicating that something was wrong with the photo- 
metric gravities (Section l23] ). 

As a second step, we calculated Fe I and Fe II abundances for 
a set of models with parameters extending more than 3cr around 
the photometric estimates of Table |3] i.e., about ±300-500 K in 
Tetf, ±0.3-0.6 dex in logg and ±0.5 km s"' in v,, depending on 
the star. This large grid of calculated abundances was used to 
refine our photometric estimate of the atmospheric parameters. 



Parameter 


Here 


TOO 


Y05 


Resolution 


30000 


30-60000 


28000 


S/N 


70-100 


>100 


30-100 


T.//(K) 


4650 


4730 


4700 


logg (dex) 


2.8 


2.4 


2.3 


V, (km s"') 


1.3 


1.8 


1.3 


[Fel/H] 


0.06±0.01(±0.10) 


-0.01±0.02 


-0.01+0.02 


[Fell/H] 


0.01±0.03(±0.10) 


-0.01±0.05 


0.01+0.03 


[ff/Fe] 


0.01±0.01(±0.07) 


0.01+0.03 


0.12+0.06 


[Al/Fe] 


0.06±0.06(±0.05) 


0.08±0.01 


0.16+0.03 


[Ba/Fe] 


0.26±0.05(±0.06) 


0.07±0.00 


0.02+0.00 


[Ca/Fe] 


-0.13+0.02(+0.03) 


0.09+0.05 


0.09+0.02 


[Co/Fe] 


0.11+0.02(±0.06) 


0.05+0.07 


0.01+0.05 


[Cr/Fe] 


0.01±0.03(±0.06) 


0.12+0.05 


— 


[La/Fe] 


0.06±0.05(±0.09) 


-0.06+0.00 


0.13+0.03 


[Mg/Fe] 


0.29±0.03(±0.10) 


0.11+0.00 


0.18+0.02 


[Na/Fe] 


0.10±0.02(±0.04) 


0.25+0.00 


0.24+0.06 


[Nd/Fe] 


0.01±0.29(±0.10) 


— 


— 


[Ni/Fe] 


0.06±0.02(±0.03) 


0.05+0.04 


0.06+0.03 


[0/Fe] 


-0.05+0.09(±0.08) 


0.04+0.00 


0.10+0.04 


[Sc/Fe] 


~0.02±0.08(±0.08) 


0.09+0.05 




[Si/Fe] 


0.09±0.02(±0.08) 


0.11+0.04 


0.11+0.04 


[Til/Fe] 


-0.07±0.02(±0.07) 


0.04+0.05 


0.05+0.03 


[Till/Fe] 


-0.07±0.08(±0.06) 


0.11+0.07 




[V/Fe] 


0.13±0.04(±0.14) 


0.09+0.03 




[Y/Fe] 


-0.04±0.02(±0.12) 


-0.12+0.08 





We chose the model that satisfied simultaneously (within the 
uncertainties^) the following conditions: (i) the abundance of 
Fe I lines should not vary with excitation potential jt',,.,; (ii) the 
abundance of Fe I lines should not vary significantly with EW, 
i.e., strong and weak lines should give the same abundanc^; 
(Hi) the abundance of Fe I lines should not differ sigificantly 
from the abundance of Fe II lines; (iv) the abundance of Fe I 
lines should not vary significantly with wavelength. 

Using the spreads of the Fe I and Fe II abundances of vari- 
ous lines, and the uncertainties in the slopes of the above condi- 
tions, we estimated the typical Icr uncertainties on the spectro- 
scopic parameters: about ±100 K on T^ff, ±0.2 dex on logg and 
±0.1 km s"' on v,. 

The resulting spectroscopic parameters (Table [3]) were al- 
ways in good agreement with the photometric ones, within the 
quoted uncertainties, with a tendency of the spectroscopic grav- 
ities to be systematically higher by 0.3-0.5 dex (see above). 
However, these differences are always easily accomodated 
within the uncertainty ranges (due to photometric errors, uncer- 
tainties in distance moduli and reddenings, and bolometric cor- 
rections). 



4.2. Abundance Calculations 

Abundance calculations and spectral synthesis (for oxygen) 
were done using the last updated version of the abundance cal- 
culation code originall y described by[S pite ( 1967). We used the 
model atmospheres by lEdvardsson et aL (J 993) . We also made 



We basically considered a slope consistent with zero when the 3 cr 
spread around the fit was larger than the maximum [Fe/H] difference 
implied by the fitted slope at the extremes of the interval covered by the 
abscissae (be it^„, EW or A). 

We decided not to use the iMag^ ( Il984h effect, because we pre- 
fer to have internally consistent abundances from each line, and the 
difference between the two methods in our case appears small (Av, < 
0.2 km s-')- 
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use of ABOMAN, a tool developed by E. Rossetti at the INAF, 
Bologna Observatory, Italy, that allows for the semi-automatic 
processing of several objects, using the above abundance cal- 
culation codes. ABOMAN performs automatically all the steps 
needed to choose the best model, and provides all the graphical 
tools to analyze the results. 

When the best model was found for each star, abundances 
and abundance ratios of all the species of interest were cal- 
culated, by averaging the results for each line of that element 
(Table|7J. Abundance ratios were always computed with respect 
to Fe I. Cluster averages were computed as weighted averages 
of the results for each star in the cluster and, if necessary, of dif- 
ferent ionization stages of each element (Table |8]l. In all Tables, 
[a/Fe] is the weighted average of all a-elements abundances. 

We can compare our results for star 141 in M 67 (see also 
Sectio n [33] ) with the ones by Tautvaisiene ( 2000) and^Yong e t alj 
(l2005h . We find a general agreement for both atmospheric pa- 
rameter and abundance ratios, with the only significant excep- 
tion of barium (but see the discussion in Section 16.31 ). calcium 
(discussed in Section |6^ and titanium. For all the discussed el- 
ements our results, as discussed in Section|6] are in broad agree- 
ment with the whole body of high resolution abundances for OC, 
so we consider our results sound. 

4.3. Internal Abundance Uncertainties 

Random uncertainties in the EW measurement process and in the 
\oggf determinations were taken into account by averaging the 
abundances determinations obtained for different lines and us- 
ing cr/ yfiniines) as the final internal error These are reported in 
Table|2] and they are of the order of ~0.01 dex for Fe I which had 
the highest number of lines, and could reach up to ~0.2-0.3 dex 
for elements relying on a handful of lines such as Nd, for exam- 
ple. Additional systematic (from line to line) and random (from 
star to star) uncertainties in EW measurements, due to contin- 
uum placement, had a negligible impact on the final abundances 
in our specific case (Section[33]l. 

For the spectral synthesis of oxygen, the uncertainty of the fit 
was computed by bracketing the best fitting synthetic spectrum 
with two spectra with altered abundance. The bracketing spec- 
tra were chosen to overlap the Icr poissonian uncertainty on the 
observed spectrum. The abundance difference of the two brack- 
eting spectra with the best fit were averaged together to produce 
the estimated uncertainty, reported in Table |7] (and in our solar 
analysis. Table [TOb. 

4.4. Uncertainties due to the Choice of Stellar Parameters 

The choice of stellar parameters implies systematic (from line 
to line) and random (from star to star) uncertainties on the final 
abundances. Usually, to estimate the impact of the stellar param- 
eter choice on the derived abundances, each parameter is altered 
by its estimated uncertainty and the resulting abundance differ- 
ences with respect to the best model parameter abundance set are 
summed in quadrature to obtain a global uncertainty. We applied 
this method to our coolest (508 in NGC 2099) and warmest star 
(2129 in Cr 110) and obtained for the various abundance ratios 
uncertainties ranging from 0.05 to 0.45 dex, with a typical value 
around ~0.10 dex. 

However, as noted bv Cavrel et al.l (|2004|) . this method pro- 
duces a very conservative estimate of the uncertainty, because it 
neglects the natural correlation among stellar parameters occur- 
ring when the so called "spectroscopic method" (Section |4. It 



Table 9. Uncertanties due to the choice of stellar parameters. 



Ratio 


NGC2099-508 
(T=4500 K) 


Crl 10-2129 
(T=4950 K) 


Average 


[Al/Fe] 


+0.04 


±0.07 


±0.05 


[Ba/Fe] 


+0.04 


±0.08 


±0.06 


[Ca/Fe] 


±0.03 


±0.03 


±0.03 


[Co/Fe] 


±0.07 


±0.05 


±0.06 


[Cr/Fe] 


±0.09 


±0.03 


±0.06 


[Fe/H] 


±0.16 


±0.04 


±0.10 


[La/Fe] 


±0.10 


±0.07 


±0.09 


[Mg/Fe] 


±0.11 


±0.09 


±0.10 


[Na/Fe] 


±0.03 


±0.06 


±0.04 


[Nd/Fe] 


±0.09 


±0.12 


±0.10 


[Ni/Fe] 


±0.02 


±0.03 


±0.03 


[0/Fe] 


±0.10 


±0.07 


±0.08 


[Sc/Fe] 


±0.07 


±0.09 


±0.08 


[Si/Fe] 


±0.12 


±0.03 


±0.07 


[Ti/Fe] 


±0.07 


±0.06 


±0.06 


[V/Fe] 


±0.21 


±0.07 


±0.14 


[Y/Fe] 


±0.15 


±0.08 


±0.12 



is adopted. Covariance terms should therefore be included to 
properly take into account such depe ndencies among the pa- 
rameters (see IMc William et al.L I1995L for a detail ed treatment 
of th e problem). The practical method proposed bv lCavrel et al.l 
1*2004) assumes that Tgff largely dominates the abundance results 
and, therefore, Tefr has to be varied by its estimated uncertainty 
(^100 K in our case. Section 1431 1. A new "second best" model 
has to be identified with the new value of Tgff by varying v, and 
\ogg accordingly, to minimize as much as possible the slopes 
of the relations described in Section 1431 This method naturally 
and properly takes into account both the main terms of the error 
budget and the appropriate covariance terms. 

We therefore altered the temperature of our hottest and 
coolest stars (see above) by -i-lOO K and by -100 K. We re- 
optimized all the parameters and re-calculated all the abundance 
ratios. The final uncertainties are the average of the uncertainties 
calculated with the higher and lower temperature and are shown 
in Table |9] The average between the uncertainties of these two 
cases is taken as a reliable estimate for the impact of the choice 
of stellar parameters on our abundance ratios. We added these 
external uncertainties between parenthesis after each abundance 
ratio and we summed them in quadrature with the internal errors 
to produce the errobars in each Figure. 

4.5. Other Sources of Uncertainty 

The following additional sources of systematic uncertainties are 
not explicitly discussed here, but should be taken into account 
when comparing our abundance estimates with other works in 
the literature: 

- systematic uncertainties due to the choice of the solar ref- 
erence abundances, which are not discussed here. Our abun- 
dances can be reported to any solar reference abundance with 
the information in Table [TOl 

- uncertainties due to the choice of log g f values, which can 
be estimated by comparing our log g f values with other lit- 
erature values (see Section [331 i: 

- uncertainties in the whole analysis due to more sophisticated 
effects such as, NLTE, HFS, isotope ratios, that are difficult 
to evaluate in some cases (these are mentioned whenever 
known or relevant in Sections [3] and|6j 
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Table 7. Abundance Ratios for sigle cluster stars, with their internal uncertainties (Section |4.3l l. For external uncertainties see 
Tables 







Cr 110 




NGC 2099 (M 37) 






NGC 2420 




Ratio 


Star 2108 


Star 2129 


Star 3144 


Star 067 


Star 148 


Star 508 


Star 041 


Star 076 


Star 174 


[Fel/H] 


0.02±0.01 


0.05+0.01 


0.02+0.01 


0.01+0.01 


-0.03+0.01 


0.07+0.01 


-0.07+0.01 


-0.06+0.01 


-0.03+0.01 


[Fell/H] 


0.00±0.06 


-0.04+0.08 


0.00+0.06 


-0.01+0.02 


-0.07+0.04 


0.05+0.06 


-0.09+0.06 


0.00+0.02 


-0.07+0.04 


[«/Fe] 


-0.02+0.02 


0.01+0.02 


0.04+0.02 


-0.04+0.01 


0.00+0.01 


0.03+0.01 


-0.03+0.01 


-0.03+0.01 


0.03+0.02 


[Al/Fe] 


-0.02±0.04 


-0.06+0.08 


-0.03+0.06 


-0.06+0.08 


-0.10+0.07 


-0.04+0.07 


-0.10+0.04 


-0.10+0.05 


-0.09+0.05 


[Ba/Fe] 


0.48±0.02 


0.49+0.04 


0.52+0.06 


0.60+0.05 


0.55+0.04 


0.58+0.04 


0.54+0.04 


0.58+0.05 


0.65+0.02 


[Ca/Fe] 


-0.08+0.04 


-0.04+0.06 


-0.09+0.07 


-0.06+0.05 


-0.08+0.04 


-0.09+0.04 


-0.08+0.03 


-0.11+0.04 


-0.08+0.05 


[Co/Fe] 


0.13+0.05 


-0.08+0.04 


0.01+0.04 


-0.04+0.02 


-0.04+0.02 


-0.05+0.02 


0.10+0.05 


-0.02+0.03 


0.05+0.03 


[Cr/Fe] 


0.00+0.04 


0.06+0.06 


0.03+0.08 


-0.01+0.06 


0.01+0.05 


0.05+0.05 


-0.02+0.03 


-0.09+0.06 


-0.02+0.05 


[La/Fe] 


0.11+0.09 


0.12+0.03 


-0.15+0.04 


0.08+0.02 


0.16+0.08 


0.13+0.08 


0.23+0.01 


0.12+0.09 


0.07+0.13 


[Mg/Fe] 


0.19+0.07 


0.01+0.14 


0.16+0.08 


0.28+0.05 


0.26+0.04 


0.27+0.03 


0.09+0.05 


0.10+0.09 


0.14+0.08 


[Na/Fe] 


-0.08+0.02 


-0.06+0.08 


0.00+0.03 


0.08+0.05 


0.10+0.04 


0.05+0.08 


-0.13+0.09 


-0.04+0.06 


0.02+0.07 


[Nd/Fe] 


0.05+0.16 


0.29+0.13 


0.51+0.28 


0.23+0.29 


0.26+0.25 


0.33+0.32 


0.12+0.16 


0.40+0.28 


A 1 O . XT A 

0.18+N.A. 


[JNl/reJ 


-0.02+0.02 


-0.01+0.02 


-0.04+0.03 


-0.04+0.03 


-0.02+0.02 


A A 1 1 A AT 

-0.01+0.02 


A AQ 1 A AT 

0.1)3+0.02 


A A 1 1 A AO 

-0.01+0.02 


A A 1 1 A AO 

-0.01+0.02 


LU/i<ej 


0.08+0.08 


-0.07+0.12 


0.02+0.14 


0.12+0.07 


0.25+0.15 


U.ZZ+U. I J 


A Al _i_A 1 f. 
— U.U1±U. 10 


A OA-i-Ct 1 A 
U.Z4+U. 14 




[Sc/Fe] 


0.11+0.05 


-0.07+0.06 


-0.14+0.10 


-0.05+0.05 


-0.01+0.07 


A AT 1 A 1 T 

-0.02+0.12 


A f\n 1 A A 1 

0.07+0.01 


A 1 1 1 A AT 

0.11+0.0/ 


0. lo+O.Oo 


rc; /c^i 
Lai/reJ 


0.02+0.03 


0.04+0.02 


0.06+0.03 


0.07+0.03 


0.1 1+0.02 


A A/^ 1 A AQ 

0.06+0.03 


A AC 1 A AT 
0.05 + 0.02 


A AC 1 A AT 

0.05+0.02 


A AT 1 A Al 

0.02+0.03 


L ili/Fej 


0.01+0.03 


0.00+0.03 


0.01+0.04 


-0.10+0.02 


-0.10+0.02 


n 07-1-0 09 


0^^-t-O 09 

— U.UDitU.UZ 


n nQ-t-O 09 
— U.U7±U.UZ 


u.uo±u.u^ 


[ ilil/rej 


0.01+0.13 


-0.04+0.07 


0.11+0.05 


-0.01+0.06 


0.00+0.05 


O 0'^-l-0 OC 


O 0/1 -uO 9T 


O 07-1-0 OQ 


u. iz+u.u / 


LV/i<eJ 


0.28+0.06 


-0.03+0.05 


-0.09+0.07 


-0.09+0.02 


-0.08+0.04 


O OT-i-O 0/1 


A 1 1 .A A^ 


n nQ-uO nc 


0. 13+0.04 


Li/reJ 


0.02+0.12 


-0.17+0.08 


-0.03+0.18 


-0.03+0.23 


0.09+0.09 


U.US+U.ZU 


A 1 A-i-r\ 1 A 


n n/^-i-O 09 
U.uD+U.uz 


u.oz+u.oy 




NGC 2682 (M 67) 




NGC 7789 










Ratio 


Star 141 


Star 223 


Star 286 


Star 5237 


Star 7840 


Star 8556 








[Fel/H] 


0.06+0.01 


0.04+0.01 


0.06+0.01 


-0.01+0.01 


0.03+0.01 


0.12+0.01 








[Fell/H] 


0.01+0.03 


0.02+0.05 


-0.05+0.04 


-0.06+0.05 


-0.03+0.07 


0.07+0.06 








[tt/Fe] 


0.01+0.01 


-0.03+0.01 


-0.05+0.01 


-0.02+0.02 


-0.02+0.02 


-0.08+0.02 








[Al/Fe] 


0.06+0.06 


0.02+0.05 


0.03+0.05 


-0.01+0.10 


0.05+0.11 


-0.13+0.10 








[Ba/Fe] 


0.26+0.05 


0.27+0.05 


0.23+0.05 


0.49+0.06 


0.36+0.07 


0.46+0.03 








[Ca/Fe] 


-0.13+0.02 


-0.20+0.03 


-0.15+0.02 


-0.10+0.04 


-0.12+0.07 


-0.24+0.03 








[Co/Fe] 


0.11+0.02 


0.05+0.04 


0.17+0.04 


-0.02+0.04 


-0.03+0.02 


-0.01+0.03 








[Cr/Fe] 


0.01+0.03 


-0.01+0.04 


0.06+0.05 


0.05+0.04 


-0.02+0.03 


0.16+0.45 








[La/Fe] 


0.06+0.05 


-0.06+0.06 


0.07+0.03 


0.14+0.06 


0.06+0.07 


0.14+0.15 








[Mg/Fe] 


0.29+0.03 


0.20+0.08 


0.23+0.05 


0.26+0.07 


0.24+0.04 


0.13+0.06 








[Na/Fe] 


0.10+0.02 


-0.02+0.04 


0.22+0.08 


0.04+0.01 


0.09+0.08 


-0.15+0.01 








[INa/rcJ 


0.01+0.29 


0.10+0.25 


0.10+0.21 


0.21+0.38 


-0.17+0.27 


U. J J±U.ZU 








[Ni/Fe] 


0.06+0.02 


0.04+0.02 


0.04+0.02 


-0.02+0.02 


0.00+0.02 


-0.02+0.03 








[0/Fe] 


-0.05+0.09 


0.09+0.11 


0.13+0.10 


0.21+0.15 


0.13+0.10 


0.23+0.20 








[Sc/Fe] 


-0.02+0.08 


-0.05+0.04 


0.01+0.05 


0.10+0.09 


0.06+0.10 


0.08+0.12 








[Si/Fe] 


0.09+0.02 


0.12+0.02 


0.06+0.02 


-0.01+0.02 


-0.02+0.03 


0.03+0.04 








[Til/Fe] 


-0.07+0.02 


-0.11+0.02 


0.00+0.02 


-0.02+0.03 


-0.03+0.03 


-0.14+0.04 








[Till/Fe] 


-0.07+0.08 


0.02+0.06 


-0.06+0.13 


-0.07+0.29 


-0.19+0.17 


0.14+0.05 








[V/Fe] 


0.13+0.04 


0.06+0.04 


0.25+0.04 


0.05+0.04 


-0.01+0.04 


-0.16+0.06 








[Y/Fe] 


-0.04+0.02 


-0.06+0.05 


0.00+0.03 


0.22+N.A. 


-0.01+N.A. 


0.07+0.18 









- small additional uncertainties due to the particular choice of 
atmospheric models and of the abundance calculation code. 



4.6. The Sun 

To test the whole abundance determination procedure, including 
EW measurement, choice of lines and atomic parameters, and 
the uncertainties determination criteria, we performed an abun- 
dance analysis of the Sun, and checked that we obtain solar val- 
ues for all elements, within the uncertainties. We used the so- 
lar spectrum from the ESO spectrograph HARPS, in La Silla, 
Chile, obtained by observing Ganymed43- The spectral reso- 
lution, R^45000, is comparable to ours, while the S/N^350 is 
much higher 

To measure EWs, we used DAOSPEC and the same linelist 
used for our programme stars. We then compared our solar 

http://www.ls.eso.org/lasilla/sciops/3p6/harps/monitoring/sun.html 



EWs with two different lit erature sets, the first by 'Moor e et al.l 
(Il966h and the second by iRutten &v an der Zalm (1984]). The 
median difference between our EW and the ones by Moore 
is EWo/io-EWMoo/-e=0.9 mA (with an inter-quartile range of 
+2. 7 mA), based on 225 lines in common, while the one with 
the iRutten & van der Zaknl (11984 ) measurements is EWoao- 
EWrz84=-0.5 mA (with an inter-quartile range of ±2.1 mA), 
based on 112 lines in common. For completeness, we note 
that EWMoo;e-EWRZ84=1.5 mA (with an inter-quartile range of 
±2.8 mA), based on 390 lines in common. We are therefore sat- 
isfied with our solar EW measurements. 

We then performed our abundance analysis (as in 
Section 14.2b by exploring the following atmospheric parame- 
ters ranges: reff,0=57OO-58OO, in steps of 50 K; logg0=4.3- 
4.5, in steps of 0.1 dex; and v,=0.5-1.5, in steps of 0.1 km s"'. 
The resulting best model has reff.0=575O K; logg0=4.4 dex; 
and Vf-0.8 km s~', i n good agreement with the accepted val- 
ues dAndersenl 1 19991) . Our adopted reference solar abundances 
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Table 8. Average cluster abundances, obtained with the weighted average of the single stars abundances. 



Ratio 


Cr 110 




NGC 2099 


NGC 2420 




M67 




NGC 7789 


[Fe/H] 
[a/Fe] 


+0.03±0.02(d 
+0.01±0.03(d 


:0.10) 
:0.07) 


+0.01±0.05(+0.10) 
0.00±0.04(±0.07) 


-0.05+0.03(±0.10) 
-0.02±0.03(±0.07) 


+0.05: 
-0.02: 


t0.02(: 

t0.03(H 


tO.lO) 
=0.07) 


+0.04±0.07(±0.10) 
~0.04±0.03(±0.07) 



[Al/Fe] 

[Ba/Fe] 

[Ca/Fe] 

[Co/Fe] 

[Cr/Fe] 

[La/Fe] 

[Mg/Fe] 

[Na/Fe] 

[Nd/Fe] 

[Ni/Fe] 

[O/Fe] 

[Sc/Fe] 

[Si/Fe] 

[Ti/Fe] 

[V/Fe] 

[Y/Fe] 



-0.04±0.02(H 
+0.49±0.02(= 
-0.07±0.02(d 
+0.01±0.10(= 
+0.02±0.03(= 
+0.03±0.18(= 
+0.16±0.07(= 
-0.06±0.05(d 
+0.23±0.20(= 
-0.02±0.01(d 
+0.03±0.08(= 
+0.01+0.13(= 
+0.04±0.02(= 
+0.02±0.04(= 
+0.05±0.19(= 
-0.10±0.12(d 



0.05) 

:0.06) 

0.03) 

:0.06) 
:0.06) 
:0.09) 
rO.lO) 

0.04) 

:0.10) 
:0.03) 
:0.08) 
:0.08) 
:0.07) 
:0.06) 
:0.14) 
:0.12) 



-0.06±0.03(H 
+0.57±0.02(= 
-0.08±0.01(d 
-0.04±0.01(d 
+0.02±0.03(= 
+0.09±0.05(= 
+0.27±0.01(= 
+0.09±0.02(= 
+0.27±0.05(= 
-0.02±0.01(d 
+0.15±0.08(= 
-0.03±0.03(d 
+0.09±0.03(= 
-0.08±0.04(d 
-0.08±0.03(d 
+0.07±0.06(= 



0.05) 

:0.06) 
:0.03) 

0.06) 

:0.06) 
:0.09) 
:0.10) 
:0.04) 
:0.10) 
:0.03) 
:0.08) 

0.08) 

:0.07) 

0.06) 
0.14) 

:0.12) 



-0.10±0.01(H 
+0.62±0.07(= 
-0.09±0.02(d 
+0.03±0.06(= 
-0.03±0.03(d 
+0.23±0.09(= 
+0.09±0.06(= 
-0.04±0.07(d 
+0.19±0.17(= 
+0.00±0.02(= 
+0.15+0.18(= 
+0.07±0.05(= 
+0.04±0.01(= 
-0.04±0.08(d 
-0.05+0. 13(d 
+0.05±0.08(= 



0.05) 

:0.06) 
:0.03) 
:0.06) 

0.06) 

:0.09) 
rO.lO) 

0.04) 

:0.10) 
:0.03) 
:0.08) 
:0.08) 
:0.07) 

0.06) 
0.14) 

:0.12) 



+0.03±0.02(= 
+0.25±0.02(= 
-0.16±0.03(d 
+0.08±0.06(= 
+0.01±0.03(= 
+0.05±0.06(= 
+0.27±0.04(= 
+0.08±0.09(= 
+0.08±0.05(= 
+0.05±0.01(= 
+0.04±0.10(= 
-0.03±0.04(d 
+0.10±0.02(= 
-0.04±0.06(d 
+0.15±0.13(= 
-0.05±0.04(d 



:0.05) 
:0.06) 
:0.03) 
:0.06) 
:0.06) 
:0.09) 
:0.10) 
:0.04) 
:0.10) 
:0.03) 
:0.08) 

0.08) 

:0.07) 

0.06) 

:0.14) 

0.12) 



-0.03±0.09(= 
+0.45±0.05(= 
-0.18±0.09(H 
-0.02±0.01(H 
+0.01±0.05(= 
+0.11±0.05(= 
+0.22±0.07(= 
-0.05±0.13(H 
+0.17±0.30(= 
-0.01 ±0.0 1(H 
+0.16±0.06(= 
+0.08±0.02(= 
-0.01±0.02(H 
-0.03±0.09(H 
-0.01±0.09(H 
+0.08±0.09(= 



:0.05) 
:0.06) 
;0.03) 
:0.06) 
:0.06) 
:0.09) 
rO.lO) 
;0.04) 
:0.10) 
;0.03) 
:0.08) 
:0.08) 
;0.07) 
;0.06) 
:0.14) 
:0.12) 



Table 10. Solar abundance values. 



Element 


[X/H]rf(,r,-, 


■eel 


logeGS96 


log6A05 


Al 


-0.34±N.A.(: 


t0.05) 


6.47+0.07 


6.37+0.07 


Ba 


+0.34±0.02(: 


t0.06) 


2.13+0.05 


2.17+0.07 


Ca 


-0.09±0.03(= 


t0.03) 


6.36+0.02 


6.31+0.04 


Co 


+0.05±0.02(: 


t0.06) 


4.92+0.04 


4.92+0.08 


Cr 


+0.02±0.05(: 


t0.06) 


5.67+0.03 


5.64+0.10 


Fel 


+0.01±0.01(: 


tO.lO) 


7.50+0.04 


7.45+0.05 


Fell 


-0.03±0.02(= 


tO.lO) 


7.50+0.04 


7.45+0.05 


La 






1.17+0.07 


1.13+0.05 


Mg 






7.58+0.05 


7.53+0.09 


Na 


-0.10±0.01(= 


t0.04) 


6.33+0.03 


6.17+0.04 


Nd 






1.42+0.06 


1.45+0.05 


Ni 


+0.01±0.01(: 


t0.03) 


6.25+0.01 


6.23+0.04 


O 


-0.01±0.07(d 


:0.08) 


8.87+0.07 


8.66+0.05 


Sc 


+0.06±0.02(: 


t0.08) 


3.17+0.10 


3.05+0.08 


Si 


-0.08±0.01(= 


t0.07) 


7.55+0.05 


7.51+0.04 


Ti 


+0.03±0.05(: 


t0.06) 


5.02+0.06 


4.90+0.06 


V 


-0.14+0.02(= 


t0.14) 


4.00+0.02 


4.00+0.02 


Y 






2.24+0.03 


2.21+0.02 



jGrevesse et al 1 Il996h are shown in Table [TOl along with the 
abundance ratios derived as described. As can be seen, all the 
derived abundance ratios are compatible with zero, within the 
uncertainties, with the exception of Al and Ba. For Al, only one 
(/l=6698A) very weak (EW=18 mA) line could be measured in 
the solar spectrum, while we used about 8 lines for the analy- 
sis of our red clump giants. The lines at 6696 amd 6698A are 
known to give tendentially lower values than the other Al lines 
(iReddv et al.Ll200llGratton"etan.l200l . so we do not worry too 
much that our lone 6698AE3 solar line gives a low [Al/Fe] result. 
Ba is discussed in Section |63] For some elements (La, Mg, Nd, 
Y) no ratio could be determined either because their lines ap- 
pear too weak in the sun, or because the solar spectrum range 
(5000-7000A) does not contain the lines we used in this paper. 



" Due to the shorter spectral range of the solar spectrum and to a 
spectral defect, we could only measure one Al line in the Sun. 



5. Cluster-by-Cluster Discussion 

5.1. Cr110 

Collinder 110 is a poorly studied, intermediate-age OC lo- 
cated at Q'j2ooo=06:38:24 and 5j2ooo=+02:01:00. We could not 
find any high resolution spectroscopic study of this cluster in 
the literature, but photom etric s t udies have been c o nducte d by 
Tsarevskii & Abakumo^ (|1971|) . iDawson & lannal (1 19981) and 
Bragagli a & Tosil 020031) . Reddening, distance and ages deter- 
mined by these authors are included in Table |4] Concerning 
metallicity, while the fir st two studies assume solar metallicity, 
'Bragagha&^s3 (l2003l) find, as a result of their synthetic dia- 
gram analysis, two equally good solutions, one at solar metal- 
licity and the other at slightly sub-solar metallicity (Z=0.008). 
A re-evaluation of the same data by .Bragaglia & Tosi (.20061) 
favours the slightly subsolar value. 

Low resolution spectroscopy using the infrared calcium 
triplet by Carrera et al. (2007) gave: [Fe/H] =-0.01 ±0.07 dex 
in the Carretta & Gratton ( 199% scale, [Fe/H]=0.0±0.3 dex 
in the iZinn & West (1984) scale and [Fe/H]=-0.19±0.21 dex 
in the [Kraft & Ivans (2003) scale. Our determination of 
[Fe/H]=-i-0.03±0.02 (±0.10) dex is in good agreement with all 
these estimates, given the large uncertainties involved in pho- 
tometric and low-resolution spectroscopic metallicity estimates. 
The other element ratios determined here have no previous liter- 
ature values to compare with, but the comparisons in Section |6] 
show that they behave like expected for a solar metallicity OC. 
Our radial velocity estimate for Cr 110, <Vr>=41.0+3.8 km s"' 
(Section I2.2[ l, is ir i good agreement with the CaT value by 
ICarrera et all (l2007l) of 45±8 km s ' . 

5.2. NGC 2099 (M 37) 

NGC 2099 (M 37) is located in the Galactic anti-center direc- 
tion in Auriga q'j2ooo=05:52:18 and (5j2ooo=+32:33:12. Since it 
is near and it appears as a relatively rich and large cluster, it 
has been photometrically studied by several authors to derive 
accurate magnitudes, proper motions, and a census of variable 
stars (for historical references see Kalirai et al., 2001). All the 
papers that derived reddening, distance and age are also cited in 
Table |4] Photometric studies generally attribute a solar metal- 
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Table 11. Literature abundance determinations for M 67 based on high resolution spectroscopy. 
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TOO" 




soo* 


Y05^ 


R06'' 


P08'' 


S09' 


R=A/6A 




30000 


30-60000 


30000 


28000 


45000 


100000 


50000 


S/N 


50-100 


>100 


40-70 


30-100 


90-180 


^80 


100-300 


[Fe/H] 


+0.05±0.02( 


tO.lO) 


-0.03+0.04 


-0.05: 


tO.Ol 


+0.02+0.01 


+0.03+0.01 


+0.03+0.04 


0.00+0.01 


[Al/Fe] 


+0.03±0.02( 


t0.05) 


+0.14+0.07 






+0.17+0.01 


-0.06+0.04 


-0.03+0.11 


— 


[Ba/Fe] 


+0.25±0.02( 


t0.06) 


+0.07+0.12 


+0.21: 


t0.20 


-0.02+0.03 


. . 


. . 





[Ca/Fe] 


-0.16±0.03(: 


t0.03) 


+0.05+0.09 


-0.02: 


t0.03 


+0.07+0.02 


+0.05+0.04 


+0.03+0.07 





[Co/Fe] 


+0.08±0.06( 


t0.06) 


+0.08+0.07 









. 








[Cr/Fe] 


+0.01+0.03( 


tO.06) 


+0.10+0.09 






— 


-0.01+0.04 


+0.03+0.09 


— 


[La/Fe] 


+0.05±0.06( 


t0.09) 


+0.13+0.10 






+0.11+0.02 











[Mg/Fe] 


+0.27±0.04( 


tO.lO) 


+0.10+0.06 


-0.11: 


t0.08 


+0.16+0.02 


-0.01+0.03 


— 


— 


[Na/Fe] 


+0.08±0.09( 


t0.04) 


+0.19+0.08 


-0.05: 


t0.05 


+0.30+0.03 


+0.04+0.07 


-0.02+0.07 




[Ni/Fe] 


+0.05±0.01( 


tO.lO) 


+0.04+0.08 


+0.04: 


t0.06 


+0.08+0.02 


-0.01+0.04 


-0.02+0.07 




[O/Fe] 


+0.04+0. 10( 


tO.03) 


+0.02+0.06 


-0.02: 


t0.08 


+0.07+0.02 


-0.04+0.01 


-0.07+0.09 




[Sc/Fe] 


-0.03±0.04(: 


t0.08) 


+0.10+0.06 














[Si/Fe] 


+0.10±0.02( 


t0.07) 


+0.10+0.05 






+0.09+0.03 


+0.03+0.04 


-0.03+0.06 




[Ti/Fe] 


-0.04±0.06(: 


t0.06) 


+0.04+0.13 






+0.12+0.02 


-0.03+0.04 


-0.02+0.11 




[V/Fe] 


+0.15±0.13( 


t0.14) 


+0.15+0.16 














[Y/Fe] 


-0.05±0.04(: 


t0.12) 


+0.01+0.14 















" iTautvaisiene! i200(f), from 6 red clump stars. 

' [Shetrone & Sandquis| i2000h . from 4 turn-off stars (we ignored blue stragglers). 

[Vong et al. ( 2005), from 3 red clump stars. 

[kandich et al. (2006), from 8 dwarfs and 2 slightly evolved stars. 
' iPace et al.l j2008l) . from 6 main sequence stars. 
^ I Santos et alJ (l2009l) . from 3 red clump giants and 6 dwarfs. 



licity to NGC 2099 (e.g., iMermilhod et all Il996h . MetalUcity 
estim ates based on photometry can only be found i n Janes et al] 



([1988 ). who give [Fe/H] =0.09 dex: [Marsha ll et alJ ([20051). who 
give [M/H]=+0.05 +0.05 and iKalirai &T osi (2004), who give 
Z<0.02. 

Surprisingly, when considering the wealth of photometric 
studies, M 37 lacks any specific low or high resolution study 
aimed at determining its chemical composition. Our values 
therefore fill this gap, and show that in all respects this cluster 
has a typical solar metallicity, with all element ratios close to 
zero within the uncertainties. On the other hand, radial velocity 
determinations for this cluster are quite abundant (Section [2. 2| i 
and appear in good agreement with our determination. 



5.3. NGC 2420 

NGC 2420 (q'j2ooo=07:38:23 and 5j2()oo=+21:34:24) has always 
been considered the definitive example of the older, mod- 
erately metal-deficient OC beyond the solar circle. Several 
good quahty imaging s tudies appea r ed already in the 60s and 
70s ('Sarmafc Walkei', '1962'; 'West*, '1967b"; 'Cannon & Llovd', 
[1970; van Altena & Jones, 1970; McClure et al., 1974, 1978, 
to name a few), and more recent photometries appeared 
in a variety of photome t ric sy stems (the most cited be- 
ing lAnthonv-Twarog et all Il990l) . Its intermediate status be- 
tween the solar-metallicity OC near the sun and the clearly 
metal-deficient population of globular clusters tagged it early 
on as a potential transition object between the two popula- 
tions, with metallicity determinations — both photometric and 
spectroscopic — placi ng it at an [Fe/HI va l ue around the 
one of 47 Tii c (e.g. Pilachowski et al., '198Ct ICohenl Il980t 
ICanterna efsli [1986; Smith & Suntzeff, 198j). More recent 
photometric works give some what higher TFe/H] values, rang- 
ing from -0.5 to -0.3 dex (e.g. iTwarog et all 119971 ; [Friel et al.L 



12002); I Anthonv-Twarog etakl l2006l) . but still significantly lower 
than the value of [Fe/H] =-0.05+0.03 (+0.10) dex we found here. 

However, both lCohenI (Il980h and iPilachowski et all (Il980h 
noted that NGC 2420 should be significantly more m etal rich 
than the Globular Clusters they analyzed, i.e., M 71 ( Cohenl 
1980) and 47 Tuc (iPilachowski et all 1 19801) . by some ^0.5 dex. 
Since they placed M 71 and 47 Tuc around [Fe/H]=-1.3, they 
consequently placed NGC 2420 at [Fe/H] =-0.6. The resolu- 
tion of their spectra (R< 10000) was much lower than ours, 
but if we trust their analysis in a relative sense, and consider 
more recent metallicity estimates for 4 7 Tuc and M 71 (-0.76 
and -0.73 respectively, iHarrisL Il996h . we would then place 
NGC 2420 around [Fe/H]:^-0.1 or -0.2. Having said that, it 
is surprising that there are no modern high resolution stud- 
ies of a cluster that was considered so important in the past. 
The highest spectral resolution empl oyed to study NGC 2420 
is R:^ 15000 fSmith & Suntzefj,ll987l) . with a spectral coverage 
of o nly 200A, giv ing [Fe/H] =-0.57. Only the preliminary work 
of Freeland et al.l ([2(302) has suggested a higher, slightly subso- 
lar [Fe/H] value for NGC 2420. We also note that our [Fe/H] 
brings NGC 2420 more in line with other OC in the Galactic 
trends discussed in Section [7] Also, we cannot ignore the simi- 
larity with the case of NGC 7789 (Sectionl5^.5l ). where high reso- 
lution spectroscopy by Tautvaisiene et alJ ([2005) and us provides 
a much higher abundance than the previous photometric and 
low/medium-resolution studies. Clearly, further high resolution 
spectroscopy with modern instruments, possibly with R^50000 
and S/N^lOO is needed for this cluster 



5.4. NGC 2682 (M 67) 

Among the old OC, M 67 (q'j2ooo=08:51:18, (5j2()oo=+1 1:48:00) 
is quite close to us, with low reddening (Table |4]i and solar 
metallicity, so it is one of the most studied ope n clusters, and 
a good target to look for solar twins and analogs (iPasquini et al.L 
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>0 



Table 12. Literature abundance determinations for NGC 7789 
based on high resolution spectroscopy. 
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[Ca/Fe] 


-0.18±0.09(: 


t0.03) 


+0.14±0.07 


[Co/Fe] 


-0.02±0.01(: 


t0.06) 
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[Cr/Fe] 


+0.01±0.05(: 


t0.06) 


-0.05±0.09 


[Mg/Fe] 


+0.22±0.07(: 


tO.lO) 


+0.18±0.07 


[Na/Fe] 


-0.05+0. 13(: 


t0.04) 


+0.28±0.07 


[Ni/Fe] 


-0.01±0.01(: 


tO.lO) 


-0.02±0.05 


[0/Fe] 


+0.16±0.06(: 


t0.03) 


-0.07±0.09 


[Sc/Fe] 


+0.08±0.02(: 


tO.08) 


-0.02+0.07 


[Si/Fe] 


-0.01±0.02(: 


tO.07) 


+0.14+0.05 


[Ti/Fe] 


-0.03+0.09(: 


t0.06) 


-0.03+0.07 


[V/Fe] 
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+0.09+0.12 


[Y/Fe] 
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" iTautvaisiene et al.l ( |2005|) . from 6 giants and 3 red clump giants. 



120081) . Since the first pioneering studies at the beginning of XX 
century, we have a few hundred papers published to date (see 



iBurstein et aLlll986l:ICarraro et al.L[l99"6l:lYadav et alll2008l for 

more references). Therefore, we have included M 67 in our sam- 
ple because it acts as a fundamental comparison object, that en- 
ables us to place our measurements in a more general frame- 
work. 

Among the vast literature on M 67, there are sev- 
eral deterr ninations o f its metallicity , with various meth- 
ods (e g. iDemarauel 119801; ICohenl [T980t iFov & Prousd. 
19811 Ijanes&S mithL [19841 lBurstein"etaLl [1984 1198^ 



Brovvnl 
19901: 



Il987t: iG arcia Lopez et al.', '1988'; 'Cavrel de Strobel' 
Hobbs & Thorburn^ 1991; Friel & Boesaaard, 1992; 
Friel & Janes, 1993; Janes & PhelDsJl994t iFriel et all l200l 
Balaguer-Nunez et al.> ,2007; .Marshall et all 12001 " to name a 



few), all typically converging to a solar value, high resolution 
abundance determinations have been derived for both giants and 
dwarfs, with many studies devoted to light elements such as 
lithium and beryllium and their implicat i ons for mixing theori es 
(iPasquini et al.Lll997t[jones et al.Lll999l;lRandich et al.[l2007h . 

Table [TT| shows a comparison of our results with 
some of the mo st rece n t high resolution (R>20000) deter- 
minations CTautvaisiend 200C ; Shetrone & Sa ndguist ._ 20001: 
iRandich et al.L 120061; iPace et al.L 120081; ISantosT t al.," 20091). The 
overall comparison is extremely satisfactory for all elements, ex- 
cept maybe for Mg, Na, Ba and Ca (see also Section [TT)- For 
Mg, Na and Ba the large spread in literature demonstrates the 
difficulties in measuring these elements. For Ca we see that our 
value is marginally lower than other literature deterrminations. 
As explained in Section [3T1 this is most probably due to the 
large uncertainties on the Calcium log gf values. 

5.5. NGC 7789 

NGC7789 (q'j20()o=23:57:24 and 5j2ooo=+56:42:30, or /=1 15.53 
and ^7^-5. 39) is a rich and intermediate-age OC, with a 
well defined giant branch, a well-populated main-sequence 
turnoff, and a substantial population of blue stragglers 
(McNamara, 1980; Twarog & Tyson, 1985; Milone & Lathamj, 
Il994l) . Several photomet ric stud ie s ha v e been carried 
out (some examples are iKustneii 119231; iReddishl 119541; 



Bur bidge & Sandagel [19581; iJanesl [19771: [Martinez Roger et all 
11 994;'Jahn et al.','1995'; 'Gim et al.', '1998a'; 'Vallenari et al.','2000^, 
[Bartasiute & Tautvaisiene, 2004; Bramich et al., 2005) and its 
parameters are reasonably well known. 

Abundance determinations through photometry and 
low/medium-resolution spectroscopy all give sub- solar values 
around [Fe/HW).2 (Pilachowski, 1985; Friel^&Janesi [19931; 
fSchonb erner et al.L 2001; Friel et al., 2002), i.e., much lower 
than our fFe/H1^0.04±0.0 7 (+0.10) dex. However , a more 
conforting comparison with iTautvaisiene et al.[ ( [2005h is shown 
in Table [T2| Their spectra have resolution and S/N similar to 
ours, and most abundance ratios in common show an excellent 
agreement. Minor discrepancies arise for some elements such as 
Ca (but see the discussions in Sections [3 . 1 I and [5r4] i, Al (but they 
used only one doublet while we used four), Na and O. Since 
they do not list their log gf values, and other ingredients of 
the abundance analysis were similar to ours, we cannot explain 
the Na-O discrepancies, but we suspect that they must be due 
log gf differences. 



6. Abundance Ratios Discussion 

We compare our abundance ratios with data from the literature, 
assembled as follows. For the Milky Way fi eld stars, we use the 
Thi ck and Thin Disc m easurements from iReddv et al.l (|2003[) 
and iReddv et al.l (l2006l) . who performed homeogeneous abun- 
dance calculations of a few hundred F/G dwarfs selected from 
the Hipparcos catalogue. We added abundance ratios, based on 
high resolution spectroscopy, for 57 old OC from various litera- 
ture sources (Table[T3]l. When more than one determination was 
available for one cluster, we simply plotted them all to give a 
realistic idea of the uncertainties involved in the compilation. 



6.1. Iron-peak Elements Ratios 

When compared with literature (Figure [4]l, our iron-peak ele- 
ments appear all solar and in good agreement with the results for 
the Disc and other OC. In particular, cobalt and chromium have 
the best agreement and smallest spreads. Although Sc, V and Co 
are known to possess HFS that may lead to an increased scat- 
ter and an overestimated ratio, they do not appear significantly 
different from solar for our target stars, so we did not attempt 
any detailed HFS analysis. Nevertheless, the effect of increased 
scatter and overestimated abundance are visible, expecially for 
vanadium, both in our data and in the Discs stars, as well as in 
the other clusters from the literature. 

A puzzling effect is seen in Figure [4| in the [Ni/Fe] 
ratio. All the data from Disc stars are very close to so- 
lar (<[Ni/Fe]>=-0.02±0.02), and so are our determinations 
(<[Ni/Fe]>=0.00±0.03), but the other OC high resolution data 
appear slightly enhanced (< [Ni/Fe] >=0. 06+0. 04), lying system- 
atically above the Disc ones. Such a ^0.05 dex offset is well 
within the uncertainties of abundance determinations in gen- 
eral, but since it appears systematic in nature, we are still left 
without a clear explanation. Our [Ni/Fe] ratios anyway are a bit 
lower than the other OC determinations, although still compati- 
ble within the uncertainties. 



6.2. a-elements Ratios 

We obtained abundances of Ca, Mg, O, Si and Ti. As can be 
seen from Figure[5] Si and Ti appear practically solar, within the 
respective uncertainties, and in very good agreement with litera- 
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Table 13. Literature sources and [Fe/H] values for high resolution (R> 15000) based abundance ratios of old OC. 



Cluster 


[Fe/H] 


Reference 


Cluster 


[Fe/H] 


Reference 


Be 17 


-0.10 


Friel et al. (2005) 


NGC 2141 


-0.26 


Yong et al. (2005) 


Be 20 


-0.61 


Yong et al. (2005) 




0.00 


Jacobson etal. (2009) 




-0.30 


Sestito et al. (2008) 


NGC 2158 


-0.03 


Jacobson et al. (2009) 


Be 22 


-0.32 


Villanova et al. (2005) 


NGC 2243 


-0.48 


Gratton & Contarini (1994) 


Be 25 


-0.20 


Carraro et al. (2007) 


NGC 2324 


-0.17 


Bragagha et al. (2008) 


Be 29 


-0.44 


CaiTaro et al. (2004) 


NGC 2360 


+0.07 


Hamdani et al. (2000) 




-0.18 


Yong et al. (2005) 




+0.04 


Smilianic et al. (2008) 




-0.31 


Sestito et al. (2008) 


NGC 2420 


-0.57 


Smith &Suntzeff (1987) 


Be 31 


-0.40 


Yona et al. (2005) 


NGC 2447 


+0.03 


Hamdani et al. (2000) 


Be 32 


-0.29 


Bragaglia et al. (2008) 


NGC 2477 


+0.07 


Bragagha et al. (2008) 


Be 66 


-0.48 


Villanova et al. (2005) 




-0.01 


Smilianic et al. (2008) 


Be73 


-0.22 


Carraro et al. (2007) 


NGC 2506 


-0.20 


Carretta et al. (2004) 


Be75 


-0.22 


Carraro et al. (2007) 


NGC 2660 


+0.04 


Bragagha et al. (2008) 


Blanco 1 


+0.04 


Ford et al. (2005^ 


NGC 3532 


+0.04 


Smilianic et al. (2008) 


Cr 261 


-0.22 


Friel et al. (2003J 


NGC 3680 


-0.04 


Pace et al. (2008) 




-0.03 


Carretta et al. (2005) 




+0.04 


Smilianic et al. (2008) 




-0.03 


De Silva et al. (2007) 


NGC 3960 


+0.02 


Bragaglia et al. (2008) 




+0.13 


Sestito et al. (2008) 


NGC 5822 


+0.04 


Smilianic et al. (2008) 


Hyades 


+0.13 


Sestito et al. (2003) 


NGC 6134 


+0.15 


Carretta et al. (2004) 




+0.13 


Paulson et al. (2003) 




+0.12 


Smilianic et al. (2008) 




+0.13 


De Silva et al. (2006) 


NGC 6253 


+0.46 


Carretta et al. (2007) 


IC 2391 


-0.03 


Randich et al. (2001) 




+0.36 


Sestito et al. (2007) 


IC 2602 


-0.05 


Randich et al. (2001) 


NGC 6281 


+0.05 


Smilianic et al. (2008) 


IC 2714 


+0.12 


Smilianic et al. (2008) 


NGC 6475 


+0.14 


Sestito et al. (2003) 


IC 4756 


-0.15 


Jacobson et al. (2007) 


NGC 6633 


+0.07 


Smilianic et al. (2008) 




+0.04 


Smilianic et al. (2008) 


NGC 6791 


+0.40 


Peterson & Green (1998) 


IC 4651 


+0.11 


Carretta et al. (20()4) 




+0.35 


Origlia et al. (2006) 




+0.10 


Pasquini et al. (2004) 




+0.39 


Carraro et al. (2006) 




+0.12 


Pace et al. (2008) 




+0.47 


Carretta et al. (2007) 


M 11 


+0.10 


Gonzalez &l Wallerstein (2000) 




+0.30 


Boesgaard et al. (2009) 


M34 


+0.07 


Schuler et al. (2003) 


NGC 6819 


+0.09 


Bragaglia et al. (2001) 


M67 


-0.03 


Tautvaisiene (2000) 


NGC 6939 


0.00 


Jacobson et al. (2007) 




-0.01 


Yong et al. (2005) 


NGC 7142 


+0.08 


Jacobson et al. (2008) 




+0.03 


Randich et al. (2006) 


NGC 7789 


-0.04 


Jautvaisiene et al. (20()5) 




+0.03 


Pace et al. (2008) 


Pleiades 


-0.03 


Randich etal. (2001) 


Mel DO 


-0.38 


Gratton & Contmni Q994) 




+0.06 


(jeuran cl Moniei (zOUo) 




-0.33 


Sestito et al. (2008) 


Praesepe 


+0.04 


Friel & Boesgaard (1992) 


Mel 71 


-0.30 


Brown etal. (1996) 




+0.12 


Pace et al, (2008) 


NGC 188 


+0.01 


Randich et al. (2003) 


Rup4 


-0.09 


Carraro et al. (2007) 


NGC 1817 


-0.07 


Jacobson et al. (2009) 


Rup7 


-0.26 


Carraro et al. (2007) 


NGC 1883 


-0.20 


Villanova et al. (2007) 


Saurer 1 


-0.38 


Carraro et al. (2004) 




-0.01 


Jacobson et al. (2009) 


Tom 2 


-0.45 


Brown et al. (1996) 


NGC 2112 


-0.10 


Brown et al. (1996) 




-0.28 


Frinchabov et al. (2008) 



ture determinations for both the Discs stars and the other OC. O, 
Ca and Mg give instead marginally discrepant enhancements. 

For O, we note that the spread both in our data and in the 
literature is greater that in any other or-element. This is partly 
due to the well known problems of determining O from the 
6300A lone line, or the 6363A weak line, or from the IR triplet 
at 7770A, that requires NLTE corrections. IVloreover, some old 
literature work uses Solar reference abundances reaching as high 
as 8.93, which can explain some of the lowest [O/Fe] literature 
estimates. Given the large spread, the tendency of our [O/Fe] 
measurements to lie on the upper envelope of the other OC high 
resolution data is probably irrelevant. Solar O enhancements 
would probably be more in line with the other a-elements, while 
the sub-solar values found generally in literature point towards 
Wolf-Rayet as additional contributor s of O, with a stronge r 
metallicity dependence of the O yields dMcWilham et al.Ll2008h . 

In the case of Ca, our values are marginally inconsistent with 
the bulk of field and OC literature determinations. A few lit- 
erature measurements of OC ratios are however as low as our 



values. These inconsistencies could be explained with the large 
uncertainties in the literature logg/ values for Calcium lines 
(~0.2 dex, see discussion in Section [TT]). Given these large ad- 
ditional uncertainties, we finally concluded that [Ca/Fe] is basi- 
cally compatible with solar values in all the clusters examined. 

Concerning Mg, we know already that the log values of 
some lines are still not very we ll determined (Section [TT])- We 
also know dGratton et all 1 19991) that some lines require NLTE 
corrections. We could find no correction factors for the lines we 
were able to meas ure in our spectra, but we noticed that those 
lines examined by Gratt on et al.l (1 19991) which hss&Xex similar to 
our lines, require NLTE corrections of about +0. 1-0.5 dex. This 
correction would make our [Mg/Fe] values even higher, reach- 
ing an enhancement of 0.2-0.6 dex with respect to solar. Another 
possibility is that our lines have a non negligible HFS, because 
they are dominated by odd isotopes, but we could find no further 
information in the literature. We could only notice that other au- 
thors find such relativ ely high values of [Mg/Fe] in OC (such as 
lBragagliaetan.l2008h . 
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Fig. 4. Comparison between our iron-peak element results (large 
black dots), the high resolution measurements of other OC listed 
in Table [13] (large dark grey do ts) and field stars belonging to 
the Thin Disc (light grey dots, iReddv et all [2003) and to the 
Thick Disc (tiny light grey dots, Re ddv et all [2006). Errorbars 
on our results are the quadratic sum of internal uncertainties and 
uncertainties due to the choice of stellar parameters (SectionHJl. 



When the average [a/Fe] values are calculated, however, all 
the programme stars and the cluster averages appear perfectly 
compatible with solar, within relatively small uncertainties, as 
expected (see Tables |7] and |8]i. The [a/Fe] ratio is also discussed 
further in Section]?] 

6.3. Heavy Elements Ratios 

We measured the heavy s-process elements Ba, La and Nd and 
the light s-process element Y. La does not require any synthesis 
to take into account HFS, since the three lines we used are al- 
ways in the linear part of the curve of growth, and in fact Figure|6] 
shows good agreement with the sparse literature values. Yittrium 
and neodimium are also in agreement with the literature data, al- 
though the measurements of Nd in OC are scarce and scattered. 
Our values of [Nd/Fe] have a tendency of being towards the up- 
per envelope of the Dies stars, but this is not significant if we 
consider the large uncertainties involved (see also the discussion 
in Section |3T2] ). 

Concerning Ba, we find high values both in our programme 
stars and in the sun itself (Se ction[43]l. The s ame result has been 
found by other authors (e.g.. lBragaglia et allbOOSl) . While a de- 
tailed study of the barium abundance is out of the scope of the 



Fig. 5. Comparison between our a-elements ratios and the liter- 
ature ones. Symbols are the same as in Figure|4] 



present paper, we note that very recently iD'Orazi et al. 1 (I2009h 
used a detailed HFS analysis of barium in OC and showed that 
the overabundance can be thus reduced roughly by ~0.2 dex. 
Looking at Figure |6] we see that in fact our [Ba/Fe] are on the 
high side of the OC data, w hich in turn have a huge spread. Data 
from lD'Orazi et all (I2009h . who revised the Ba abundances for 
20 OC using spectral synthesis to take HFS into account, are to- 
wards the lower envelope of the OC abundances (open stars in 
Figure |6]). As can be seen, some [Ba/Fe] ehnancement remains 
in their high quality data, that is apparently well correlated with 
the cluster ages (see their Figure 1). Still no clear explanation 
is available, since the current evolution models and yields do 
not reproduce the data correctly at young ages, where the en- 
hancement is higher and more uncertain (up to [Ba/Fe]^0.6 dex 
for ages around 10*^). In summary, most of the [Ba/Fe] enhance- 
ment we see in our measurements should be due to HFS effects, 
but some part of it c ould be real (up to 0.2 dex, see Fig. 2 by 
D'Orazi et al.. |2009|) . A hint of a descending slope of [Ba/Fe] 
in OC appears , that is not appaeren t among Disc stars. Further 
studies such as lD'Orazi et al.ld2009l) are necessary on large sam- 
ples of Disc and OC stars to get firmer constraint on the chemical 
evolution of Ba in the Galactic Dies. 

6.4. Ratios of Na and Al and Anticorrelations 

We also derived Na and Al abundances, since these elements are 
quite easy to measure in OC and there is a vast body of liter- 
ature measurements to compare with. Figure |7] shows that for 
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Fig. 6. Comparison between our s-process elements ratios and 
the literature ones. Symbols are the same as in FigureH] except 
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Fig. 8. A search for (anti)-correlations of Al, Mg, Na and O 
among our targets stars. The four panels show different planes 
of abundance ratios, where stars belonging to each cluster are 
marked with different symbols. Dotted lines show solar values, 
solid lines show linear regressions and the typical uncertainty 
(~0. 1 dex) is marked on the lower right corner of each panel. 
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Al there is a general agreement between our data, the Disc ra- 
tios and the OC ratios, although there is a large scatter in the 
literature data. NLTE corrections for the four Al doublets and 
the ty pe of stars studied he r e cou ld not be found in the liter- 
ature. iBaumtiller & GehrenI (Il997h give corrections for hotter 
(Teff>5000) and higher gravity (logg>3.5) stars, that suggest 
that either the corrections are negligible, or they are slightly neg- 
ative at lower temperatures and gravities. 

In the case of Na, the spread in the literature data is even 
larger, and there are both data points with significant Na en- 
hancement and points with typical solar values. Part of the 
scatter depends on the need for NLTE corrections. According 
to lGratton etall (199^ . the 5682-5688A and the 6154-6160A 



doublets at EW^lOO mA both require corrections of about 
<0. 05-0. 10 dex, for solar stars like the ones considered here. 
The NLTE corrected abundances should be higher than the LTE 
uncorrected ones: this should make our [Na/Fe] LTE measure- 
ments in better agreement with literature measurements. If the 
observed enhancement in [Na/Fe] should prove to be real for 
OC, this would set O C stars completely apart form Disc stars 
dDe Silva et al.L l2009l) . If the large spread will also prove to be 
intrinsic, this would suggest the possibility of light elements 
chemical anomalies similar, although much less pronounced, to 
the ones found in Globular Clusters. 

In fact Al and Na, together with Mg and O (and C and N) 
show puzzling (anti)-con'el ations in almost all of the studied 
Galactic Globular Clusters iGratton et al.L |2004 . The most in- 
teresting fact concerning the chemical anomalies, is that they 
have never been found outside globular clusters. They are not 
present in the field populations of the Milky Way and its sur- 
rounding dwarf galaxies, and they only recently have been found 
in Forna x and LMC clusters ( Letarte et al.L 120061: iJohnson et all 
2006; M ucciarelh et all l2009l) . Rece nt dedicated searches fo r 
anti-correlations in OC are the one bv'M artell & SmithI ( l2009h . 
on the strength of CH and CN bands in NGC 188, NGC 2158 and 
NGC 7789, the one bv lSmilianic et al.1 (|2008|) . based on high res- 
olution sp ectroscopy of C, N, O, Na and '^C/'^C and the one by 
iDe Silva et al. (2009), that compiles and homogenizes literature 
Na-O high resolution data. No clear-cut sign of anti-correlation 
has been found yet. 

If we build the usual (anti)-correlation plots for our five OC 
(Figure [8]), we find no clear sign of chemical anomalies, and in 
all four plots the spread of each ratio is still compatible with 
the typical uncertainty of our measurements, which is of the or- 
der of 0.1-0.2 dex, depending on the element. In particular, in 
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the [0/Fe]-[Na/Fe] plane, all 15 star s roughly occupy the solar 
region around zero that in Figure 5 by Carretta et al. (2006) con- 
tains normal stars only (see also De Silva et al., 2009). A possi- 
ble exception to this total absence of correlations is the [Na/Fe]- 
[Al/Fe] plane, where a hint of a correlation can be noticed. 
Statistically, this is not significant and small variations in Teff 
could induce a similar weak correlation. For this reason, anti- 
correlations are usually a more robust si gn of chemical a noma- 
lies. Still, in light of the discussion by ISmilianic et al.l (l2008l) 
about the Na-0 anti-correlation, our Na-Al results is suggestive. 
If further studies will show that some chemical anomalies of 
these elements are present in Galactic OC, our analysis (together 
with that of S milianic et al.L i2008) shows that they must be of a 
much smaller extent than in Globular Clusters: 0.2-0.3 dex at 
most in the [Al/Fe], [Mg/Fe], [Na/Fe] and [O/Fe] for the kind 
of clusters studied here. In any case, the lack of relations in OC 
would point towards one or more of the following environmental 
causes for the presence of anticorrelations in globular clusters: 
(i) relatively low metallicity (below solar); (ii) dense environ- 
ment; (Hi) total cluster mass of the order of -IO'^Mq or more; 
(iv) undisturbe d environment (e.g. away from the Disc, see also 
ICarrettall2()06b . 



0.5 



-0.5 



-0.5 




7. Galactic Trends 

As said earlier, OC are the fundamental test particles for the 
study of the chemical evolution of the Galactic Disc, and as 
such, they produce two of the strongest constraints on chemi- 
cal models: the Galactic radial trends and the Age-Metallicity 
Relation (AMR). Since a careful homogeneization of litera- 
ture data (including not only element ratios, but also ages and 
Galactocentric radii) is out of the scope of the present paper, 
we have used the literature data of Table [T3j averaging all esti- 
m ates for a single cl uster together. We complemented with data 
bv lFriel et al.l (l2002h for those OC lacking high resolution mea- 
surements. Using 28 OC in common between the two datasets, 
we found an average offset of [ Fe/H] =-0. 1 6 ±0. 1 3 dex, in the 
sense that the measurements by iFriel et all (|2002|) are on av- 
erage smaller tha n the o nes based on high resolution. We cor- 
rected iFriel et al] (|2002|) data by this amount befo re plotting 
them . We extracted OC ages from the compilation of iDias et al.l 
(l200Zfy^ and for the Ga lactocentric Radii (Rgc) we used pri- 
marily Friel et al"] (|2002|) . complemented by the WEBDA, and 
filled in the few missing clusters with data from the papers of 
Table [13] The resulting radial trends and AMR for [Fe/H] and 
[a/Fe] (computed as in Section |4!2] | are plotted in Figures |9] and 
[TOland discussed below. 



7.1. Trends with Galactocentric Radius 

The trend of abundances with Galactocentric radius Rgc gives 
a strong constraint for the models of Galactic chemical evo- 
lution as far as the Disc formation mechanism is concemaJB 



It is now widely accepted (see for example Qw arog et al 
fT997t IFriel etsll l200l lYong et al.L 12001 ISestito et al 



2008; 



We are aw are that at least in the case of NGC 6791, the age given 
by iDias et all {2002) is quite different from other literature estimates 
(citations in Tablsll3l>. being lower by at least 2 Gyr. However, building 
a homogeneous age scale is a non-trivial task, and it is beyond the scope 
of the present paper. 

" A far stronger constraint would be the variation of this trend with 
age. Such a study is at the moment not possible, given the small number 
of clusters studied with high resolution in a homogenous way. 



Fig. 9. Trends of [Fe/H] (top panel) and [a/Fe] (bottom panel) 
with Gal actocentric Radius . Light grey dots in the top panel are 
OC from IFriel et all (|2002|) . grey dots are the OC compiled in 
Table IT3] and black dots are o ur data. NGC 2420 estimates by us 
and Smi th & Suntzefil (1 19871) are connected with a thick dotted 
line. A linear fits with uncertainty is drawn across all points in 
both panels. 



iMagrini et al.L l2009l) that there is a clear trend of decreasing 
metallicity, measured as [Fe/H], with increasing Roc- Such a 
trend is clearly detected not only in OC, but also in field stars 
(B-stars and Cepheids), H II regions and planetar y nebulaj^ 
(see Chia ppini et al., 2001 : lAndrievskv et~al .. 2004; Yong et aTl 



2005; L emasle et all 120081 for some review of literature data). 

The first large studies of homogeneous OC abundances 
(summarized in the review byEiiiri995), found that old OC 
(older than the Hyades) were extending outwards in the Disc, 
much farther out than young OC, and they found a well defined 
slope out to Rgc- 16 kpc. The spread around this slope was at 
the time around ~0.2 dex, i.e., generally compatible with (or per- 
haps slightly larger than) the measurement uncertainties. Such a 
slo pe comes na t urally in most Galactic chemical evolution mod- 
els (iTbsi L IT9821 Il988l: iMatteucci & FrancoisL IT989I: iTbsiL [1996 : 
Chiappin ietaU 120011: jAndrievs kv et al.L 120041: IColavitti et al. , 
2009t iMagrini et all 12009), when different star formation and 



infall rates are assumed for the inner and outer Disk. To repro- 
duce most of the observational constraints, a differential Disc 
formation mechanism is often assumed, either with the inner 
Disk formed first and then growing in radius (inside-out forma- 
tion) or with the whole Disk evolving simultaneously, but with 
a much more intense (and sometimes fast) evolution in the cen- 
tral, denser parts. A prediction of all models is that the metallic- 
ity gradient should change with time (although different models 
predict very different time changes, Tosi, 1996) and some predict 
that it should flatten out at larg e radii. Indee d, the first stu dies of 
anticenter and distant clusters (ICarraro et al. . 2004; Yong et al.L 



Although we know of at leas one dataset in w hich PNe show flat 
trends of oxygen and neon abundances with Rqc dStanghellini et all 
120061) . 
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120051: ICarraro et all 120071; ISestito et al.Ll2008l) showed that after 
Rgc-12 kpc, the relation flattens around a value of [Fe/H]^- 
0.3 dex. 

Not much can be said with the present data about the slope 
variation with time, but the exact value of the slope has been mat- 
ter of some debate. As said, earlier studies found a value around 
-0.09±0.01 dex kpc or-0.07 ±0.01 dex kpc ' w ith the stricdy 
homogenous measurements by Friel et al ] (| 2002|). An alterna- 
tive interpretation of a different data compilation dTwarog et al.L 
Il997i) . describes the trend as two disjoint plateaux, one around 
solar metallicity and comprising OC within Rgc -10 kpc, 
and a second one at an [Fe/H]^-0.3 outside the solar cir- 
cle. More rec ent work bas e d on high resolution compilations 
of OC data (ISestito et all l2008h find a steeper slope of - 
0.17+0.01 dex kpc ' within Rgc <14 kpc, which stifl holds 
when considering only the 10 clusters analyzed homogeneously 
by that group. A sort of bimodal distribution is observed in their 
Figure 9, where between the very steep slope of the inner clus- 
ters and the plateau of the outer ones there is a small gap almost 
devoid of OC (9<Rgc<12 kpc). 

Our results are plotted in Figure |9] where we consider the 
trend of [Fe/H] and [a/Fe] with Rgc- We do see a distinct 
slope in the inner clusters and a flattening out for the outer 
ones. However, our sampl e contains ~15 more OC than the 
one of Sestit o et alj (|2008|) . and most of them (including our 
five determinations) fall in the gap around 9<Rgc<12 kpc, 
discussed above. With the addition of these clusters we find 
a gentler slope of -0.05 ±0.01 dex kpc ', in good agree- 
ment with previous works jFriell 119951: iFriele t al., 20 0l and 
with the Disk Cepheids within Rgc-H kps (' Andrievskv et al.L 
[2004; Lemasle et al., 2008). If we exclude the clusters outside 
Rgc= 12 kpc and remove the low resolution OC from iFriel et alj 
(I2OO2I) . the slope does not steepen significantly, becoming - 
0.06±0.02 dex kpc Also, the flattening out does not seem so 
abrupt as in Figure 9 by Sestito et al. (2008). The paucity of OC 
in the flat part of the relation (we have only Be 20, Be 22, Be 29 
and Saurer 1 in our compilation) surely calls for more high reso- 
lution studies, since as the data stand now, they look compatible 
with both a plateau and a gradual change in slope. We would like 
to note that NGC 2420 (alread y discussed in S ection 15.31 ) was 
placed at [Fe/H] =-0.57 dex bv , Smith & Suntzeff (1987) , based 
on R^ 16000 spectra, while we find -0.05, in much better agree- 
ment with the global Galactic trend. This goes in the direction of 
filling the gap in the Sestito et al. (2008) compilation, and also in 
the lTwarog et alj (11997i) data-set, pointing more towards a gentle 
and continuous decrease of [Fe/H]. 

The trend of a-enhancement with Rgc is also of some impor- 
tance, since it unveils the role of SNe type la and II and their rel- 
ative contributions. Yong et al. (2005) fou nd a tendency of the a - 
enhancement to increase with Rgc, as did lMagrini et alj (l2009l) . 
who found this tendency in good ag reement with their chem ical 
evolution model. Also model A by IChiappini et al.l (1200 lb pre- 
dicted an increase of [c/Fe] with Rgc- In our compilation, the 
trend appears as a weak slope, that is still perfectly compatible 
with a flat distribution at the 1 cr level. This, together with the 
study of slope changes with time, is one typical case in which 
a high quality, homogeneous analysis of ~100 OC could give a 
clear and definitive answer. 

7.2. Trends with Age 

In spite of the fact that all models predict an evolution of Disk 
metallicity with time, albeit maybe only in the first Gyrs, and the 
fact that such a variation is observed in Disk stars dReddy et al.L 
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Fig. 10. Trends of [Fe/H] (top panel) and [a/Fe] (bottom panel) 
with cluster ages. Symbols are the same as in Figure|9] 

'2003'; 'B ensbv etal.Ll2004 . there appears to be no correlation 
at all between old OC abundances and ages (see the review by 
lFrielLll995h . More recent results did not change this picture sub- 
stantially. If confirmed, the lack of an AMR in OC would point 
towards a d i fferen t source of chemical enrichment for OC stars 
(lYong et al.L l2005l) with respect to the Disc stars. In substance, 
the metallicity of OC stars seems to be more determined by the 
location in which they formed, than by the time at which they 
formed. 

What we find here is quite encouraging, although still not 
statistically significant. We recall that our compilation includes 
57 high resolution abundance determinations, plus a handful of 
low resolution determinations by Friel et al. (2002). Although 
we have made no attempt to homogeneize the data, except for a 
-0. 16 dex correction to the low resolution abundances, this sam- 
ple i s ~50% large r than a ny compilation presented before (see 
e.g., ISestito et al.L I2008L iDe Silva et all l2009t iMagrini et all 
2OO9I) and shows that the community is proceeding fast in filling 
up the gaps of our knowledge of OC. Figure [TO] shows indeed a 
weak trend of decreasing [Fe/H] and increasing [c/Fe] with in- 
creasing age. The slopes are very gentle at best and they are still 
compatible with no trends at all. Nevertheless, if such trends ex- 
ist, we can put some constraints on them: for [Fe/H], the gradient 
should not be significantly larger than -2.6+1.1 10 " dexGyr ' 
and for [a/Fe] no larger than 1.1 ±5.0 10"" dex Gyr"'. 

8. Summary and Conclusions 

We have analyzed high resolution spectra of three red clump gi- 
ants in five OC, three of them lacking any previous high reso- 
lution based chemical analysis. Given the paucity of literature 
data, such a small sample is enough to increase the whole body 
of high resolution data for OC by ^5%. To compare our re- 
sults with the literature, we have compiled chemical abundances 
based on high resolution data of 57 clusters from the litera- 
ture. Given the recent and fast progress in the field, this sam- 
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pie is ~50% larger than previous literature compilations (e.g., 
iFriel et all llool ISestito et all 120081: iMagrini et all l2009h . The 
main results drawn by the analysis of our five clusters are: 

- We provide the first high resolution based abundance analy- 
sis of Cr 110 ([Fe/H]=+0.03+0.02 (+0.10) dex), NGC 2099 
([Fe/H]=+0.01±0.05 (±0.10) dex) andNGC 2420 ([Fe/H]=- 
0.05 ±0.03 (±0.10) dex), which only had lo w resolution de- 
terminations and the 16000 analysis by Smit h & SuntzeffI 
(1987); our new determination of the metallicity of 
NGC 2420 puts this cluster in much better agreement with 
the global Galactic trends; 

- The abundances found for NGC 7789 ([Fe/H]= +0.04±0.07 
(±0.10) dex) and M 67 ([Fe/H]= +0.05±0.02 (±0.10) dex) 
are in good agreement with past high resolution studies; 

- We provide the first high resolution based radial velocity de- 
termination for Cr 110(< y,. >=41.0±3.8kms"'); 

- We found that all our abundance ratios, with few exceptions 
generally explained with technical details of the analysis pro- 
cedure, are near-solar, as is typical for OC with similar prop- 
erties; we found solar ratios also for Na, that is generally 
found overabundant, and for O, which is generally found un- 
derabundant; 

- We do not find any significant sign of anti-correlation (or 
correlation) among Na, Al, Mg and O, in general agreement 
with past and recent results, and we can say that if such cor- 
relations indeed are present in OC, they must be much less 
extended than in Globular Clusters, amounting to no more 
than 0.2-0.3 dex at most; 

With our compilation of hterature data, we also could ex- 
amine global Galactic trends, that are extremely useful to con- 
struct chemical evolution models for the Galaxy in general and 
the Galactic Thin Disc in particular. For the metallicity gradient 
we found a slope of -0.06±0.02 dex kpc ' considering only the 
high resolution data within /?gc=12 kpc. Our compilation con- 
tains data, including our own determinations, that fill the small 
gap around 9<Rgc <12 and point towards a gentl e and continu- 
ous de crease, rather than a two s tep drop such as in|Twarog et alj 
(Il997h or a steep slope such as inl Sestito etall (|2008|) . We do find 
a flattening at Rcc > 12 kpc and a hint of an increasing [a/Fe] to- 
wards the outer Disk. Concerning the AMR, we do not find any 
strong evidence for it, and we just note some very mild trends. If 
an AMR is indeed present among OC, it must be very weak and 
we provide upper limits to its slope both in [Fe/H] and [a/Fe]. 
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